40-0182  771  FLOW  ONO  FOILU 
GEOHATERIALSIU 
OF  APPLIED  NEC 
UNCLASSIFIED  AFOSR-TR-87-88 

RE  OF  ROCKS 
CALIFORNIA 
PANICS  A 

90  AFOSR-86- 

CONCRETE  AND  OTHER 

UN IV  SAN  DIEGO  LA  JOLLA  DEP 
S  NENAT-NASSER  NOV  86 

8835  F/G  8/7 

r  1/2  T 

NL 

■ 

* 

AD-A182  771 


November  1986 


flPA/&7  7/ 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  FAGS 


1a  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


.  OSCLASSIFICATION/OOWNG RACING  SCHEDULE 


4.  FSRFORMING  ORGANIZATION  REPORT  NUM8SRIS) 


13.  OISTRIBUTI ON/AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited 


8.  MONITORING  ORGANIZATION  REPORT  NUMESR(S) 

,  AfOSR-TR-  8  7  -  U  890 

Sa.  NAME  OF  PERFORMING  ORGANIZATION  tab.  OFFICE  SYMBOL  I  7*.  NAME  OF  MONITORING  ORGANIZATION 


Department  of  AMES,  University  Wm »UeabU> 

of  California,  San  Dieeo 

Sc.  AOORBSS  (City,  State  and  ZIP  Coda) 

La  Jolla,  California  92093 


Air  Force  Office  of  Scientific  Research/NA 

7b.  AOORBSS  (City,  State  and  ZIP  Coda) 

Bolling  Air  Force  Base,  D.C.  20332 


Sa.  NAME  OF  FUNOING/SPONSORING  8b.  OFFICE  SYMBOL 

ORGANIZATION  (If  applicable! 

_  fHZdSiiltffr 

Sc.  AOORBSS  (City.  State  and  ZIP  Coda)  f 

jj-osr/  Yk\  4P 
Belling  AFB,  DC  80339 

O^r'0"'  "Flow  and  Failure  of" 

\c\Vi  *3^  (Unclassified) 

nainwnwi 

S.  Nemat-Nasser 


8b.  OFFICE  SYMBOL  8.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

ftf£dS&] Nfr  Grant  AFOSR  86-0035 _ 

10.  SOURCE  OF  FUNDING  NOS. _ 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


2302 /fS 


TASK 

WORK  UNIT 

NO. 

NO. 

C3 

13a.  TYPE  OF  REPORT 

Final  Report 

18.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED  14.  OATS  OF  REPORT  (Yr.,  Mo.,  Day)  IB.  PAGE  COUNT 

from  10^1-85  to  9-30-86  November  1986 _ 191 _ 


CQEATI  COOES _ 

GROUP  I  SUE.  OR. 


1 19.  AWT 9 ACT  (Con  tin  u*  on  rmw  if  * 


_  18.  SUBJECT  TERMS  (Continue  on  meree  If  necemery  and  Identify  by  Mock  number)  Micro— 

_ mechanics  of  flow  and  failure  of  rocks,  concrete,  and  other 

_ geomaterials;  macroscopic  nonlinear  constitutive  models; 

___  macroscopic  nonlinear  and  anisotropic  reaconsgs. 

t  and  Idem  Hfy  by  block  number) 


The  research  completed  under  this  AFOSR  award  was  the  continuation  of  efforts 
already  initiated  under  Grant  AFOSR-84-0004  at  Northwestern  University.  The  research 
aims  at  experimental  and  theoretical  investigation  of  the  micromechanics  of  flow  and 
failure  of  rocks,  concrete,  and  other  related  geo -materials  at  moderate  to  very  high 
pressures  and  temperatures.  As  a  result  of  this  work,  certain  macroscopic  nonlinear 
constitutive  models  have  been  developed,  which  reflect  realistically  the 
micromechanical  events  that  produce  observed  macroscopic  nonlinear  and  anisotropic 
responses  of  materials  of  this  kind. 

The  theoretical  work  followed  the  principal  investigator's  recent  effort  in  the 
micromechanical  modeling  of  nonlinear  material  response.  It  included  calculations  of 

(continued  on  reverse 


2G  OISTRI 8UTION/AVAI LABILITY  OF  ABSTRACT 

UN CLA88I FI 80/UNLIMITSO  O  SAME  AS  RFT.  ffl  OTIC  USERS  E] 


33a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Spencer  T.  Wu 


DO  FORM  1473, 83  APR  EDITION  OF  1  JAN  73  IS  OBSOLETE. 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22b  TELEPHONE  NUMBER 
(Include  Area  Code) 

(202)  767-4935 


22c.  OFFICE  SYMBOL 
AFOSR/NA 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAOE 


Unclassif ied 


SlCUfllTY  CLASSIFICATION  Of  THIS  FAOt 


Section  19.  ABSTRACT  (Continued) 

microcrack  initiation  under  overall  compression,  interaction  between  cracks, 
development  of  plastic  zones  and  their  interaction  with  cracks,  and  the  final  failure 
mode  of,  say,  rocks.  In  particular,  attention  was  focused  on  the  influence  of  the 
pressure  and  temperature  on  the  failure  mode  and  on  the  transition  from  brittle  to 
ductile  response.  The Modeling  of  observed  axial  splitting  and  faulting  at  moderate 
pressures  and  low  temperatures  was  completed,  and  through  some  theoretical  and  model 
experiments,  a  basic  understanding  has  been  gained  for  the  phenomenon  of  brittle - 
ductile  transition  at  elevated  pressures. 

The  experimental  effort  was  carefully  coordinated  with  the  theoretical  one.  It 
consisted  of  two  parts:  (1)  qualitative  model  studies  in  order  to  identify  and 
understand  the  involved  micromechanics;  and  (2)  quantitative  model  tests.  .'! ■ 

.  i  ' 


FINAL  SCIENTIFIC  REPORT 


for  the  period 

I  October  1985  to  30  September  1986 
on 

FLOW  AND  FAILURE  OF  ROCKS,  CONCRETE  AND  OTHER  GEOMATERIALS 
Grant  No.  AFOSR-86-0035 


S.  Neoat-Nasser ,  Principal  Investigator 
University  of  California,  San  Diego,  La  Jolla,  California  92093 


Table  of  Contents 

ga.ge 

RESEARCH  OBJECTIVES  iii 

ORGANIZATION  OF  THIS  REPORT  iv 

CHAPTER  I  1 

Research  Accomplishments  3 

1.  Introduction  3 

2.  Progress  4 

3.  List  of  Publications  Completed  10 

4.  Professional  Personnel  Associated  with  the  Research 

Project;  Degrees  Awarded  11 

5.  Interactions  (Coupling  Activities)  12 

CHAPTER  II  13 

"A  Unified  Analysis  of  Various  Problems  Relating  to  15 

Circular  Holes  with  Edge  Cracks" 

CHAPTER  III  65 

"On  Mechanics  of  Crack  Growth  and  Its  Effects  on  the  67 

Overall  Response  of  Brittle  Porous  Solids" 

CHAPTER  IV  121 

"A  Microcrack  Model  of  Dilatancy  in  Brittle  Materials"  123 

November  1986 
i 


RESEARCH  OBJECTIVES 


The  research  completed  under  this  AFOSR  award  was  the  continuation  of 
efforts  already  initiated  under  Grant  AFOSR- 84-0004  at  Northwestern 
University.  The  research  aims  at  experimental  and  theoretical  investigation 
of  the  micromechanics  of  flow  and  failure  of  rocks,  concrete,  and  other 
related  geo-materials  at  moderate  to  very  high  pressures  and  temperatures. 
As  a  result  of  this  work,  certain  macroscopic  nonlinear  constitutive  models 
have  been  developed,  which  reflect  realistically  the  micromechanical  events 
that  produce  observed  macroscopic  nonlinear  and  anisotropic  responses  of 
materials  of  this  kind. 

The  theoretical  work  followed  the  principal  investigator's  recent 
effort  in  the  micromechanical  modeling  of  nonlinear  material  response.  It 
included  calculations  of  microcrack  initiation  under  overall  compression, 
interaction  between  cracks,  development  of  plastic  zones  and  their 
interaction  with  cracks,  and  the  final  failure  mode  of,  say,  rocks.  In 
particular,  attention  was  focused  on  the  influence  of  the  pressure  and 
temperature  on  the  failure  mode  and  on  the  transition  from  brittle  to 
ductile  response.  Under  the  AFOSR-84-0004 ,  the  modeling  of  observed  axial 
splitting  and  faulting  at  moderate  pressures  and  low  temperatures  was 
completed,  and  through  some  theoretical  and  model  experiments,  a  basic 
understanding  has  been  gained  for  the  phenomenon  of  brittle -ductile 
transition  at  elevated  pressures. 

The  experimental  effort  was  carefully  coordinated  with  the  theoretical 
one.  It  consisted  of  two  parts:  (1)  qualitative  model  studies  in  order  to 
identify  and  understand  the  involved  micromechanics;  and  (2)  quantitative 


model  tests. 


ORGANIZATION  OF  THIS  REPORT 


This  report  is  organized  in  the  following  manner. 

In  Chapter  I,  we  list  the  papers  that  have  been  completed  and  provide  a 
brief  abstract  of  each  paper.  At  the  end  of  Chapter  I,  a  list  of  scientific 
articles  completed  under  this  project  is  given.  In  Chapter  I,  we  also  give 
the  list  of  participants  and  related  information.  Chapters  II  -  IV  each 
presents  a  complete  research  effort  which,  although  finished,  has  not  yet 
been  published.  Since  Grant  AFOSR-86-0035  actually  constitutes  the  final 
year  support  for  a  three -year  research  effort  which  was  started  at 
Northwestern  University  under  Grant  AFOSR-84-0004,  this  report  also  includes 
a  discussion  of  work  completed  under  Grant  AFOSR-84-0004. 
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CHAPTER  I 


RESEARCH  ACCOMPLISHMENTS 


1.  INTRODUCTION 

The  description  of  the  mechanical  response  of  rocks,  concrete,  and 
other  related  geo -materials  under  moderate  to  very  high  pressures  and 
temperatures  is  of  fundamental  importance  to  current  and  future 
technological  efforts.  Many  large-scale  computer  codes  require  accurate  and 
realistic  constitutive  modeling  of  materials  of  this  kind.  At  the  same 
time,  to  be  useful  and  effective,  these  models  must  be  simple.  Therefore, 
it  is  of  paramount  importance  first  to  develop  an  in-depth  understanding  of 
the  essential  features  responsible  for  the  behavior  of  this  kind  of  material 
in  various  loading  and  temperature  regimes,  and  then  to  seek  to  formulate 
macroscopic  constitutive  models  that  embody  the  most  essential 
micromechanical  features.  In  this  manner,  one  will  have  a  good  chance  of 
arriving  at  realistic  models  of  required  simplicity. 

One  main  feature  of  materials  at  focus  here  is  that  they  naturally 
contain  many  microdefects  and  inhomogeneities.  These  defects  and 
inhomogeneities  serve  as  stress  and  strain  concentrators,  resulting  in  local 
stresses  which  are,  in  general,  truly  three-dimensional  and  very  large,  even 
when  the  overall  applied  load  is  one -dimensional  (uni -axial)  and  relatively 
small.  The  overall  behavior,  therefore,  is  greatly  influenced  by  the 
microstructure,  and  changes  as  the  microstructure  is  changed  in  response  to 
the  applied  loads. 
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The  pressure  affects  the  cl.  '.nges  of  the  microstructure  to  a  great 
extent.  At  low  confining  pressures,  the  existing  cracks  may  open  under 
tension  or  may  undergo  frictional  sliding  when  under  overall  compression, 
with  shear  stress  acting  across  the  plane  of  the  crack.  In  either  case, 
tension  cracks  develop  at  the  tips  of  the  pre-existing  flaws,  and  grow  in 
the  direction  parallel  to  the  maximum  compression.  At  higher  pressures,  the 
growth  of  such  tension  cracks  is  soon  arrested,  but  under  overall  applied 
loads,  smaller  flaws  and  inhomogeneities  that  are  closely  spaced  begin  to 
interact.  This  leads  to  a  different  microscopic  response  and,  therefore, 
macroscopic  failure.  When  the  pressures  are  very  high,  local  plastic  flow 
is  the  major  micromechanics  of  deformation,  and  the  overall  response,  though 
highly  nonlinear  and  still  pressure-sensitive,  becomes  plastic. 

The  temperature  has  a  similar  effect.  At  low  temperatures  the  response 
is  more  brittle,  whereas  at  high  temperatures  it  is  more  ductile. 

Therefore,  macroscopic  constitutive  modeling  must  take  into  account 
such  varied  microscopic  behavior  patterns.  The  aim  of  the  research  completed 
under  this  grant  has  been  to  develop  theoretical  micromechanical  models,  to 
verify  these  by  quantitative  and  qualitative  model  experiments,  and  then  to 
seek  to  develop  appropriate  macroscopic  constitutive  relations  which  embody 
the  essential  microfeatures. 

2 .  PROGRESS 

Under  Grant  AFOSR-84-0004 ,  effort  has  been  concentrated  on  completing 
the  microscopic  modeling:  (1)  analyzing  the  phenomenon  of  faulting  under 
moderate  confining  presures;  and  (2)  understanding  the  transition  from 
brittle  to  ductile  failure,  under  suitably  high  confining  pressures,  of 
brittle  solids  containing  microdefects  such  as  pre-existing  cracks.  This 


work  was  continued  under  Grant  AFOSR-86-0035  at  the  University  of 
California,  San  Diego.  Microcracking  from  voids  in  an  elastic  solid  was 
formulated  and  solved.  In  addition,  nonlinear  constitutive  relations 
associated  with  microcracking  and  damage  have  been  developed.  In  the 
sequel,  papers  completed  and  published  under  both  grants,  A.FOSR- 84-0004  and 
AFOSR-86-0035,  are  listed  together  with  their  abstracts.  In  addition, 
wherever  appropriate,  comments  are  made  to  point  out  the  significance  of  the 
obtained  results . 


2.1  Papers  Published  Under  Grant  AF0SR-84-0004: 

2.1.1  H.  Horii  and  S.  Nemat-Nasser ,  "Elastic  Fields  of  Interacting  Inhomo¬ 
geneities,"  Int.  J.  Solids  Structures,  21  (1985)  731-745. 

In  this  paper  a  rather  effective  method  has  been  developed,  by  which 
the  interaction  between  two  or  among  several  defects  in  an  elastic 
solid  can  be  computed  on  the  basis  of  the  solution  of  a  single  defect 
in  an  infinitely  extended  homogeneous  solid.  The  results  are  funda¬ 
mental  and  allow  us  to  estimate  the  failure  mechanism  in  the  presence 
of  large  confining  pressures  where  the  interaction  among  neighboring 
flaws  leads  to  faulting.  The  abstract  of  this  paper  is  as  follows. 

A  rather  general  technique  --  called  the  "method  of  pseudotractions" 
is  presented  for  the  calculation  of  the  stress  and  strain  fields 
in  a  linearly  elastic  homogeneous  solid  which  contains  any  number  of 
defects  of  arbitrary  shape.  The  method  is  introduced  and  illustrated 
in  terms  of  the  problems  of  elastic  solids  containing  two  or  several 
circular  holes  and  solids  containing  two  or  several  cracks,  including 
the  cases  of  rows  of  holes  or  cracks.  It  is  shown  that  the  solution 
of  these  and  similar  problems  can  be  obtained  to  any  desired  degree 
of  accuracy.  Furthermore,  if  only  estimates  are  needed,  then  the 
method  is  capable  of  yeilding  closed- form  analytic  expressions  for 
many  interesting  cases,  e.g.  the  stress  intensity  factors  at  the 
crack  tips. 

2.1.2  H.  Horii  and  S.  Nemat-Nasser,  "Brittle  Failure  in  Compression: 

Splitting,  Faulting,  and  Brittle-Ductile  Transition,"  Philosophical 
Transactions  of  the  Royal  Society  of  London  319,  No.  1549  (1986) 

337-374. 

This  is  a  comprehensive  paper  which  we  hope  will  have  a  lasting 
effect  on  the  understanding  of  the  mechanics  of  brittle  failure  in 
compression  and  the  mechanism  of  brittle-ductile  transition  at 
elevated  pressures.  Considerable  effort  has  been  required  in 
preparing  and  publishing  this  paper. 


The  micromechanics  of  brittle  failure  in  compression  and  the 
transition  from  brittle  to  ductile  failure,  observed  under  increasing 
confining  pressures,  are  examined  in  the  light  of  existing 
experimental  results  and  model  studies.  First,  the  micromechanics  of 
axial  splitting  and  faulting  is  briefly  reviewed,  certain 
mathematical  models  recently  developed  for  analysing  these  failure 
modes  are  outlined,  and  some  new,  simple  closed-form  analytic 
solutions  of  crack  growth  in  compression  and  some  new  quantitative 
model  experimental  results  are  presented.  Then,  a  simple  two- 
dimensional  mathematical  model  is  proposed  for  the  analysis  of  the 
brittle-ductile  transition  process,  the  corresponding  elasticity 
boundary -value  problem  is  formulated  in  terms  of  singular  integral 
equations,  the  solution  method  is  given,  and  numerical  results  are 
obtained  and  their  physical  implications  are  discussed.  In  addition, 
a  simple  closed- form  analytic  solution  is  presented  and,  by  comparing 
its  results  with  those  of  the  exact  formulation,  it  is  shown  that  the 
analytic  estimates  are  reasonably  accurate  in  the  range  of  the 

brittle  response  of  the  material.  Finally,  the  results  of  some 

laboratory  model  experiments  are  reported  in  an  effort  to  support  the 
mathematical  models. 

In  addition,  the  following  papers  have  been  published: 

2.1.3  S.  Nemat-Nasser  and  H.  Horii,  "Micro-mechanics  of  Fracture  and 
Failure  of  Geo-materials  in  Compression,"  Proceedings  of  the  ICF  6th 
Int.  Conf.  on  Fracture,  New  Delhi,  India,  December  4-10,  1984; 
Pergamon  Press,  1984,  Vol.  1-R,  pp.  515-524. 

Recent  analytical  results  on  non-coplanar  crack  growth  in  elastic 
solids  under  far- field  compressive  stresses,  are  used  to  examine  the 
micromechanics  of  brittle  failure  in  compression.  The  three  distinct 
failure  modes  --  axial  splitting,  faulting,  and  the  transition  from 

brittle  to  ductile  response  --  observed  under  axial  compression  for 

different  confining  pressures,  are  discussed  in  terms  of  simple 
plausible  micro-mechanical  models.  The  failure  strength  and  the 
orientation  of  failure  planes,  as  well  as  the  stress  ratio  which 
marks  the  birttle-ductile  transition,  are  estimated  and  compared  with 
published  data  on  various  rocks,  arriving  at  good  correlations.  In 
addition,  certain  model  experiments  which  seem  to  support  the 
analytical  models,  are  examined. 

2.1.4  S.  Nemat-Nasser  and  H.  Horii,  "Rock  Failure  in  Compression,"  Int.  J. 
Engng .  Sci.,  Letters  in  Appl .  and  Engrg.  Sciences,  1984,  Vol.  22,  No. 
8-10,  pp.  999-1011. 

Mechanisms  of  rock  failure  in  compression  --  axial  splitting  and 
shear  failure  --  are  studied,  based  on  a  microscopic  consideration. 
Analytical  models  are  constructed  and  model  experiments  on  plates  of 
a  brittle  polymer  are  performed.  It  is  shown  that  unstable  growth  of 
tension  cracks  which  propagate  from  the  tips  of  pre-existing  cracks 
and  curve  towards  the  maximum  compressive  direction,  is  the 
fundamental  mechanism  that  produces  axial  splitting  of  a  uniaxially 
compressed  rock  specimen,  whereas  shear  failure  of  a  triaxially 
compressed  specimen  is  a  result  of  sudden  growth  of  tension  cracks  at 
tips  of  a  suitably  arranged  interacting  set  of  microcracks.  The 


simultaneous  out-of-plane  unstable  growth  of  a  suitably  oriented  row 
of  cracks  is  analyzed  and,  on  the  basis  of  this  model,  the  variations 
of  the  "ultimate  strength"  and  the  orientation  of  the  overall  fault 
plane  with  the  confining  pressure  are  estimated.  The  brittle -ductile 
transition  is  discussed  with  the  aid  of  a  model  which  includes  both 
tension  crack  extension  and  plastic  zone  development  from  the  pre¬ 
existing  cracks. 

2.1.5  H.  Horii  and  S.  Nemat-Nasser,  "Compression- Induced  Microcrack  Growth 
in  Brittle  Solids:  Axial  Splitting  and  Shear  Failure,"  J.  Geophys . 
Res.,  1985,  Vol.  90,  No.  B4,  pp.  3105-3125. 

Micromechanisms  of  rock  failure  (axial  splitting  and  shear  failure) 
are  examined  in  light  of  simple  mathematical  models  motivated  by 
microscopic  observations.  The  elasticity  boundary  value  problem 
associated  with  cracks  growing  from  the  tips  of  a  model  flaw  is 
solved.  It  is  shown  that  under  axial  compression,  tension  cracks 
nucleate  at  the  tips  of  the  preexisting  model  flaw,  grow  with 
increasing  compression,  and  become  parallel  to  the  direction  of  the 
maximum  far-field  compression.  When  a  lateral  compression  also 
exists,  the  crack  growth  is  stable  and  stops  at  some  finite  crack 
length.  With  a  small  lateral  tension,  on  the  other  hand,  the  crack 
growth  becomes  unstable  after  a  certain  crack  length  is  attained. 
This  is  considered  to  be  the  fundamental  mechanism  of  axial  splitting 
observed  in  uniaxially  compressed  rock  specimens.  To  model  the 
mechanism  of  shear  failure,  a  row  of  suitably  oriented  model  flaws  is 
considered  and  the  elasticity  boundary  value  problem  associated  with 
the  out-of -plane  crack  growth  from  the  tips  of  the  flaws  is  solved. 
It  is  shown  that  for  a  certain  overall  orientation  of  the  flaws  the 
growth  of  the  out-of -plane  cracks  may  become  unstable,  leading  to 
possible  macroscopic  faulting.  On  the  basis  of  this  model  the 
variations  of  the  "ultimate  strength"  and  the  orientation  of  the 
overall  fault  plane  with  confining  pressure  are  estimated,  and  the 
results  are  compared  with  published  experimental  data.  In  addition, 
the  results  of  a  set  of  model  experiments  on  plates  of  Columbia  resin 
CR39  containing  preexisting  flaws  are  reported.  These  experiments 
are  specifically  designed  in  order  to  show  the  effect  of  confining 
pressure  on  the  crack  growth  regime.  The  experiments  seem  to  support 
qualitatively  the  analytical  results. 

2.1.6  S.  Nemat-Nasser  and  H.  Horii,  "Mechanics  of  Brittle  Failure  in 
Compression,"  Computers  &  Structures,  1985,  Vol.  20,  No.  1-3,  pp. 
235-237. 

The  mechanics  of  brittle  failure  in  compression  is  reviewed  in  light 
of  experimental  observations  and  some  recent  results  from 
micromechanical  modeling.  The  axial  splitting,  faulting  and  the 
transition  from  brittle  to  ductile  response  are  discussed.  It  is 
concluded  that  the  descending  portion  of  the  stress-strain  curve 
(i.e.  the  strain-softening),  often  observed  for  materials  of  this 
kind,  does  not  represent  a  continuum  response:  therefore,  many 
numerical  calculations  in  the  strain- softening  range  should  be  viewed 
with  caution. 


2.2  Work  Completed  Under  Grant  AFOSR-86-0035 


Three  major  papers  and  a  short  note  were  completed,  submitted  for 
publication,  and  since  have  been  accepted.  Two  of  these  papers, 
items  2.2.1  and  2.2.2,  relate  to  the  failure  of  porous  solids  in 
compression.  The  results  apply  to  porous  rocks  and  to  ice  where 
bubbles  and  cavities  are  responsible  for  initiating  microcracking 
under  compressive  loads.  The  third  paper,  item  2.2.3,  is  a  major 
contribution  in  obtaining  overall  nonlinear  constitutive  relations 
for  brittle  solids  in  compression.  The  effects  of  microflaws  are 
formulated  and  the  corresponding  overall  response  is  computed.  The 
short  note  is  an  interesting  illustration  of  the  complex  overall 
material  response  in  the  presence  of  microcracks. 


2.2.1  "A  Unified  Analysis  of  Various  Problems  Relating  to  Circular  Holes 
with  Edge  Cracks"  Engineering  Fracture  Mechanics ,  in  press. 

A  unified  method  of  analysis  is  developed  for  various  problems 
relating  to  elastic  plates  containing  circular  holes  with  edge 
cracks.  The  method  is  based  on  the  analysis  of  a  unit  rectangular 
region  containing  a  circular  hole  with  edge  cracks,  where  the 
boundary  conditions  of  the  outer  edges  are  suitably  adjusted  in  order 
to  treat  various  problems  including  periodic  arrays  of  holes  with 
edge  cracks.  The  method  is  applied  to  five  problems,  and  accurate 
values  of  the  stress  intensity  factors  are  obtained.  These 
analytical  values  of  the  stress  intensity  factors  are  fitted  by 
polynomials  which  are  convenient  for  practical  applications. 

2.2.2  "On  Mechanics  of  Crack  Growth  and  Its  Effects  on  the  Overall  Response 
of  Brittle  Porous  Solids,"  Acta  Metallurgies,  in  press. 

This  paper  is  concerned  with  crack  growth  in  brittle  porous  solids 
under  compression  and  its  effects  on  the  overall  response  of  the 
material.  As  a  mathematical  model,  we  consider  an  elastic  solid 
containing  a  zig-zag  array  of  circular  holes  with  a  pair  of  edge 
cracks  (two-dimensional  problem),  and  solve  this  problem  by  using  a 
theory  which  gives  numerical  results  as  accurate  as  desired.  Based 
on  the  analytical  results,  we  discuss  the  crack  growth  process  and 
estimate  the  effective  Young's  moduli  as  well  as  the  stress-strain 
relation  for  porous  solids.  Our  computations  show  that  the  cracks 
emanating  from  the  poles  of  the  circular  holes  extend  in  the  axial 
direction  and  grow  --  in  most  cases  in  a  stable  manner,  but  for 
certain  cases  in  an  unstable  manner  during  an  intermediate  loading 
state  --  as  the  overall  applied  uniaxial  compression  increases, 
reaching  a  certain  limiting  maximum  length.  This  maximum  crack 
length  strongly  depends  on  the  ratio  of  the  hole  radius  to  the  hole 
spacing  in  the  loading  direction.  The  effective  Young’s  modulus  in 
the  direction  of  the  crack  growth  is  basically  determined  by  the 
initial  porosity,  and  is  little  affected  by  the  crack  length  or  its 
growth  regime,  i.e.,  whether  stable  or  unstable.  We  find  that  the 
overall  axial  stress-strain  curve  remains  monotonic,  exhibiting  no 
peak  stress  or  strain  softening,  as  cracks  extend  in  the  axial 
direction  and  reach  their  limiting  length  with  increasing  axial 
stress . 
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2.2.3  "A  Microcrack  Model  of  Dilatancy  in  Brittle  Materials,"  Journal  of 
Applied  Mechanics ,  in  press. 

For  a  solid  containing  preexisting  flaws,  overall  nonlinear 
constitutive  relations  are  developed  on  the  basis  of  a  model  which 
endows  a  preexisting  flaw  with  frictional  and  cohesive  resistance, 
and  which  includes  nucleation  and  growth  of  tension  cracks  at  the 
preexisting  flaw,  as  it  deforms  under  the  action  of  an  overall 
compressive  load.  The  preexisting  flaws  may  be  randomly  distributed 
or  may  have  an  initial  preferential  distribution.  They  may  be  of 
varying  sizes  and  orientations.  Even  when  the  flaws  are  randomly 
distributed,  their  preferential  activation,  and  the  nucleation  and 
growth  of  tension  cracks  at  preferential  flaws,  render  the  overall 
response  of  the  solid  highly  anisotropic.  As  a  first  step  toward  a 
more  complete  constitutive  micromechanical  modeling,  a  dilute 
distribution  of  preexisting  flaws  is  assumed,  rate  constitutive 
relations  are  developed  for  loading  and  unloading,  which  include 
hysteresis,  dilatancy,  and  other  characteristics  observed 
experimentally  in  rocks,  ceramics,  concrete,  and  similar  brittle 
materials.  A  number  of  illustrative  examples  are  worked  out,  and 
the  results  are  compared  to  relevant  experimental  observations. 

The  following  papers  were  revised  as  required  during  the  review 
process,  considerably  improved  and  expanded,  and  the  final  manuscripts  were 
submitted  for  publication: 


B.  Rowshandel  and  S.  Nemat-Nasser,  "Finite  Strain  Rock  Plasticity: 
Stress  Triaxiality ,  Pressure,  and  Temperature  Effects,"  Soil 
Dynamics  and  Earthquake  Engineering,  in  press. 

B.  Rowshandel  and  S.  Nemat-Nasser,  "A  Mechanical  Model  for 
Deformation  and  Earthquake  on  Strike-Slip  Faults,"  Pure  and 
Applied  Geophysics  124,  No.  3  (1986)  532-566. 


LIST  OF  PUBLICATIONS  COMPLETED 


H.  Horii  and  S.  Nemat-Nasser ,  "Elastic  Fields  of  Inhomogeneities,"  Int . 
J.  Solids  Structures,  21  (1985)  731-745. 

H.  Horii  and  S.  Nemat-Nasser,  "Brittle  Failure  in  Compression: 
Splitting,  Faulting,  and  Brittle-Ductile  Transition,"  Philosophical 
Transactions  of  the  Royal  Society  of  London  319,  No.  1549  (1986) 

337-374.  (Considerable  effort  has  been  devoted  to  completing  this 
manuscript,  which  is  a  rather  comprehensive  account  of  our  work 
supported  by  AFOSR) . 

S.  Nemat-Nasser  and  H.  Horii,  "Micromechanics  of  Fracture  and  Failure 
of  Geomaterials  in  Compression,"  Proceedings  of  the  ICF  6th  Int. 
Conf.  on  Fracture,  New  Delhi,  India,  December  4-10,  1984:  Pergamon 
Press,  1984,  Vol.  1-R,  pp.  515-524. 

S.  Nemat-Nasser  and  H.  Horii,  "Rock  Failure  in  Compression,"  Int.  J. 
Engrg .  Sci. ,  Letters  in  Appl .  and  Engrg.  Sciences ,  1984,  Vol.  22, 
No.  8-10,  pp.  999-1011. 

H.  Horii  and  S.  Nemat-Nasser,  "Compress ion- Induced  Micro-Crack  Growth 
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Res. ,  1985,  Vol.  90,  No.  B4,  pp.  3105-3125. 
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M.  Isida  and  S.  Nemat-Nasser,  "A  Unified  Analysis  of  Various  Problems 
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Metallurgica,  in  press. 
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CHAPTER  II 


A  UNIFIED  ANALYSIS  OF  VARIOUS  PROBLEMS  RELATING  TO 
CIRCULAR  HOLES  WITH  EDGE  CRACKS 

M.  Isida*  and  S.  Nemat-Nasser 

Department  of  Applied  Mechanics  and  Engineering  Sciences 
University  of  California,  San  Diego,  La  Jolla,  CA  92093 


Abstract 

A  unified  method  of  analysis  is  developed  for  various  problems 
relating  to  elastic  plates  containing  circular  holes  with  edge  cracks.  The 
method  is  based  on  the  analysis  of  a  unit  rectangular  region  containing  a 
circular  hole  with  edge  cracks,  where  the  boundary  conditions  of  the  outer 
edges  are  suitably  adjusted  in  order  to  treat  various  problems  including 
j  periodic  arrays  of  holes  with  edge  cracks.  The  method  is  applied  to  five 

I 

problems,  and  accurate  values  of  the  stress  intensity  factors  are  obtained. 
These  analytical  values  of  the  stress  intensity  factors  are  fitted  by 

polynomials  which  are  convenient  for  practical  applications. 

I 

t 

* 
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I  *  Present  address:  Department  of  Mechaical  Engineering,  Kurume  Institute 

[  of  Technology,  2228  Mukuno  Karaitsumachi ,  Kurume  830,  Japan. 
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NOMENCLATURE 


x.  y 

z 

2b,  2c 

7 

a 

a 

l 

L 

fi 

E 


Cartesian  coordinates 
complex  coordinate  -  x  +  iy 
length  and  width  of  rectangular  region 
aspect  ratio  of  rectangular  region  -  - 
hole  radius 

normalized  hole  radius  -  - 

c 

crack  length 

2 

normalized  crack  length  -  - 

a 

crack  length  parameter  -  (1+L)  ^ 
Young's  modulus  of  material 


Poisson's  ratio  of  material 
3  -  i/ 


K  - 


J"+  (plane  stress) ,  -  3  -  4i/  (plane  strain) 


shear  modulus  - 


2(l+«/) 


u,  v 


4>{z) ,  ^(z) 

^1(z),  ^(z) 

P_(s  ) 

m  m 

q_(s  ) 

m  m 


stress  components  in  Cartesian  coordinates 

reference  stress 

principal  stresses  at  infinity 

a  2 

ratio  of  principal  stresses  at  infinity  -  —  (o^><0) 

resultant  force  components  in  Cartesian  coordinates 

displacement  components  in  Cartesian  coordinates 
complex  stress  potentials 

complex  stress  potentials  in  Laurent  expansion 
densities  of  distributed  force  doublets 
weighting  functions  of  force  doublet  densities 


subdivision  number  of  crack 
subdivision  number  of  side  EF 
subdivision  number  of  side  FG 
length  of  intervals  on  side  EF 
length  of  intervals  on  side  FG 
Mode  I  stress  intensity  factor 
dimensionless  stress  intensity  factor 

dimensionless  stress  intensity  factor 
critical  stress  intensity  factor 


1 .  INTRODUCTION 


Since  the  pioneering  work  of  Griffith  [1],  it  has  bean  commonly 
accepted  that  brittle  failure  is  often  caused  by  the  formation  and  growth  of 
tension  cracks  which  usually  nucleate  at  micro- defects  such  as  preexisting 
cracks,  voids,  and  other  micro-heterogeneities.  Recent  analytical  and 
experimental  work  [2-11]  has  clearly  shown  that,  even  under  overall 
compressive  farfield  loads,  brittle  failure  is  usually  a  result  of 
nucleation  and  growth  of  micro- tension  cracks.  An  important  class  of 
micro-defects  which  may  be  a  source  of  crack  nucleation,  is  micro-cavities 
which  are  commonly  present  in  most  solids  such  as  rocks,  ceramics,  metals, 
bones,  etc.  To  be  able  to  examine  the  physical  consequences  of  strength 
degradation  by  crack  nucleation  and  growth  at  preexisting  cavities,  one 
requires  a  solution  of  the  problem  of  crack  growth  from  the  walls  of 
cavities  under  prescribed  farfield  stresses.  In  the  two-dimensional  case, 
the  elasticity  problem  of  cracks  extending  from  the  edges  of  a  circular  hole 
has  been  solved  by  Bowie  [12]  for  an  infinitely  extended  medium.  The  more 
general  problem  of  cracks  emanating  from  an  elliptical  hole  in  a  plane  has 
been  treated  by  several  authors  [13-16]  where  some  of  the  results  have  been 
put  in  a  form  convenient  for  design  purposes  [17,18], 

In  this  paper  we  consider  a  class  of  elasticity  problems  relating  to 
circular  holes  with  edge-cracks,  and  solve  these  problems  by  a  unified 
method.  This  method  is  such  that  the  results  can  be  rendered  as  accurate  as 
desired.  The  basic  problem  is  formulated  for  a  rectangular  region 
containing  a  circular  hole  with  two  edge-cracks  emanating  from  the  hole. 
Various  cases  are  then  studied  by  adjusting  the  boundary  conditions  on  the 
outer  boundaries  of  the  rectangular  region.  In  this  manner,  solutions  are 


V 

& 

•Xi 

*e* 


obtained  for  a  finite  region,  for  a  strip  containing  an  array  of  holes  with 
edge-cracks,  and  for  a  sheet  containing  doubly-periodic  holes  and  edge- 
cracks.  These  results  are  fitted  by  polynomials  which  are  convenient  for 
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practical  applications. 
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2.  THEORETICAL  ANALYSIS 

2.1  Description  of  the  problem 

Consider  a  rectangular  region  containing  a  central  circular  hole  and  a 
pair  of  equal  cracks,  emanating  from  the  hole  in  the  x-direction;  see 
Fig. 1-1  for  notation  and  definition  of  various  geometric  quantities.  A 
variety  of  problems  including  periodic  arrays  of  holes  with  edge  cracks  can 
be  analysed  by  adjusting  the  boundary  conditions  of  this  rectangular  region. 
We  consider  a  class  of  problems  for  which  the  boundary  conditions  are 
symmetric  about  both  coordinate  axes.  We  solve  this  class  of  problems  using 
suitably  symmetric  stress  potentials  defined  on  the  first  quadrant,  OEFGO  in 
Fig. 1-1,  and  by  enforcing  the  desired  boundary  conditions. 

This  paper  deals  with  the  following  five  problems: 

(1)  Uniform  normal  loads  along  EF  and  along  FG,  (Fig. 1-1). 

(2)  Infinite  strip  in  the  x-direction  with  a  periodic  array  of  holes 

and  edge  cracks  aligned  in  the  x-direction.  This  case  corresponds  to 
a  rectangular  specimen  subjected  to  uniform  displacement  u  without 

shear  along  EF,  (Fig. 1-2). 

(3)  Infinite  strip  in  the  y-direction  with  a  periodic  array  of  holes 

and  edge  cracks  aligned  in  the  x-direction.  This  case  corresponds  to 
a  rectangular  specimen  subjected  to  uniform  displacement  v  without 

shear  along  FG,  (Fig. 1-3). 

(4)  Infinite  solid  with  doubly  periodic  array  of  holes  and  edge  cracks, 

(Fig. 1-4). 

(5)  Infinite  solid  with  zig-zag  periodic  array  of  holes  and  edge 


cracks ,  (Fig. 1- 5) . 


The  conditions  on  the  outer  edges  depend  on  the  problem,  and  are 
summarized  in  Table  1.  In  this  table,  u  and  v  are  the  components  of  the 
displacement,  and  are  those  of  the  resultant  force  transmitted  across 
an  arbitrary  path  from  a  fixed  point  (Xq,  y q)  to  a  typical  point  (x,  y)  , 
and  a ^  and  a ^  are  the  applied  stresses  in  the  x-  and  y-directions , 
respectively.  The  cracks  are  assumed  to  be  traction- free  for  all  the 
treated  problems.  Hence,  the  results  are  valid  for  combinations  of  a and 
a 2  that  satisfy  this  assumption. 

2.2.  Basic  relations 


In  plane  problems  of  elasticity,  components  of  the  stress,  the 
displacement  and  the  resultant  force  with  respect  to  the  Cartesian 
coordinates  are  expressed  in  terms  of  two  complex  potentials,  $(z)  and  V’(z) . 
and  their  derivatives,  as  follows: 
a  +  a 

J  2  *  -  2  Re[*'(z)] 


+ 


ir 

xy 


z  4>"  (z)  +  0"  (z) 


(1) 


2G  (  u  -  iv)  -  k  4>(  z)-zi£'(z)-0'(z)  (2) 

P  +  iP ’  -  -  ^(z)  -  z<t>'  (z)  -  4>'  (z)  (3) 

y  x 


where  G  is  the  shear  modulus,  and  x  is  defined  in  terms  of  Poisson's  ratio, 
i/,  as  follows: 


x  -  j-  —  —  (plane  stress)  ,  -  3  -  kv  (plane  strain) 


(4) 


The  analysis  is  performed  by  superposing  two  stress  states.  The  first 


one  is  introduced  in  order  to  satisfy  the  outer  boundary  conditions,  and  is 


represented  by  the  complex  potentials  <^(z)  and  having  singularities 
within  the  hole.  The  second  stress  state  is  properly  singular  at  crack 
tips,  and  is  realized  by  distributing  suitable  force  doublets  along  the 
cracks . 

The  above  two  stress  states  are  established  in  such  a  manner  that  they 
automatically  satisfy  the  traction- free  conditions  at  the  hole.  These 
stress  states  include  unknown  parameters  which  are  then  determined  from  the 
boundary  conditions  on  the  cracks  and  the  outer  edges  of  the  rectangular 
region. 

2.3.  Displacement  and  resultant  force  due  to  <^(z)  and  V^(z) 

The  complex  potentials  z)  and  ^(z)  for  the  first  stress  state  are 
expressed  as  the  following  Laurent  series: 


^1(z) 


VlF,  z-<2"+1>+M,  z2"+1) 

n-0  2n  2n 


z-v^, 
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*{<z)  -  S[D 

1  n-0  2n 
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The  coefficients  F„  ,  M-  ,  D„  and  K„  (n  -  0,1,2, - ,  N-l)  are  real  due  to 

zn  zn  zn  zn 

the  symmetry  of  the  problem.  N  is  to  be  fixed  at  the  stage  of  the  boundary 
collocation  procedure,  as  discussed  in  Section  2.5. 

The  traction- free  conditions  at  the  hole  require  that 
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°2n- 

-<2n-1>K2n-2a"n-(2n>ZM2na4n+2 

(1  <  n  <  N-l) 

(6) 

F2n- 

•K2na4n+2-(2n+3)M2n+2  ^ 

(0  <  n  <  N- 1) 

Substituting  Eqns.  (6)  into  Eqns .  (5),  the  complex  potentials  for  the 
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first  stress  state  are  obtained  as  follows: 


,  .  .  „  .  2n+l  . .  4n  -(2n-l)  v  4n+2  -(2n+l) 

^(z)  -  Mqz  +  ZQ[M2n{z  -  (2n+l)a  z  l-K^a  z 


if/  \  2  -1  .  i  2n+l  ..  - .  4n+4  -(2n+3) 

tt( z)  -  -2Mn  a2z  +  Z  [K0  { z  -(2n+l)a  z  } 

1  u  n-0 

u  \2  4n+2  -(2n+l), 

-M2n(2n)3a  z  ] 
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The  boundary  conditions  of  the  present  problems  have  been  given  in 
Table  1  in  terms  of  stress  and  displacement  components.  But  in  the 


a 


analysis,  the  stress  conditions  are  replaced  by  those  of  the  resultant  force 
in  order  to  increase  the  accuracy  of  the  numerical  results  [19,20], 


Using  Eqns. (7)  in  Eqns. (2)  and  (3),  we  have  the  following  expression^ 


I 

I 


for  the  displacement  and  the  resultant  force  components: 

-1  N-l  9n+1  ?n 

2G(u-iv)/1,  -  M.{2a2z  *  +  (k-1)z}+  Z  M.  [kz  -(2n+l)zz 
(i)  0  n-l  2n 

,/0  .n  4n  — -  (2n-l)  ..  -2n,  ..  ,2  4n+2  -(2n+l), 

+(2n+l)a  { -kz  -(2n-l)zz  }  +(2n)2  a  z  ] 


N-l  r  2n+l  4n+2  -(2n+l)  1V_  -(2n+2), 

+  Z  K_  [-z  -a  [kz  '+(2n+l)zz  ) 

n-0  2n 


+(2n+l)a4n+4z'(2n+3) 
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y  x  (1)  0  n_^  2n‘ 
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*  V  K,  [.221>*l+a'‘n*2,2-(2..n).(2„+1),z-<2n*2)1 
n-0  ^n 

+  <2n+l)a4n+V<2r,+3)l 


(8.2) 


The  above  expressions  include  2N  unknown  coefficients  M^n  and 
(n- 0,1,2, —  ,N-1)  to  be  determined  from  the  boundary  conditions. 

2.4.  Displacement  and  resultant  force  due  to  force  doublets 


The  second  stress  state  is  realized  by  distributing  force  doublets  of 
arbitrary  density  on  the  cracks.  The  force  doublets  are  established  from 
stress  states  for  concentrated  forces  acting  in  an  infinite  solid  which 
contains  a  circular  hole,  and,  therefore,  the  traction-free  conditions  at 
the  hole  edge  are  satisfied  automatically.  Expressions  for  the  displacement 
and  the  resultant  force  components  due  to  distributed  force  doublets  have 
been  derived  by  Isida  and  Noguchi  [21],  and  are  as  follows: 
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where  Pm(sm)  (n-1,2)  are  unknown  density  functions  of  the  force  doublets, 


and  the  subscript  m  takes  1  for  the  right  and  2  for  the  left  crack, 
respectively;  the  following  notation  is  also  used  (see  Fig.  2): 


W1  “  Z  x0,m  * 

cx 

w0 - - 

2  z 

X.  w,  -  xn 

0,m  3  0,m 

X0,m 

i 

-  a  +  2 

(m-1) 

(10) 

x0,m 

i 

-  -a  -  - 

(m— 2) 

The  density  functions  Pm(sm)  in  Eqns.  (9)  are  expressed  as  [22] 

»«<*»>  -  W<  <ll> 

where  the  weighting  functions  q  (s  )  are  regular. 

mm 

Eqn.  (11)  actually  defines  force  doublet  densities  with  singularities 
at  both  ends,  s^-li/2.  It  applies  to  the  present  problem,  provided  the 
resultant  force  conditions  (instead  of  stress  conditions)  are  used,  as  can 
be  seen  by  the  accurate  results  obtained  even  with  small  number  of  of 
subdivisions  shown  in  Tables  2  to  4  of  Section  3.2. 

It  is  obvious  from  symmetry  that 

q2(-s2)  -  91(s1)  (12) 

Thus  q^(s^)  is  the  only  unknown  function  in  the  second  stress  state. 

2.5.  Boundary  conditions  based  on  resultant  force  and  displacement 

A  piecewise  linear  approximation  will  be  used  to  estimate  the  unknown 
weighting  function  q^s^).  For  the  right  crack,  this  is  shown  in  Fig.  3, 
where  equally  spaced  intervals  are  used.  Eqns.  (9.1)  and  (9.2)  now  reduce 
to  linear  expressions  in  fj  (j  -  1,2,--- ,M)  which  are  the  values  of  the 
weighting  function  at  the  ends  of  these  intervals. 


We  now  have  2N  unknown  coefficients  in  the  first  stress  state  and  M 
unknown  weights,  f  ,  in  the  second  stress  state  to  be  determined  from  the 
boundary  conditions.  Since  both  stress  states  satisfy  the  traction- free 
conditions  on  the  hole,  we  have  only  to  consider  the  boundary  conditions  on 
the  crack  and  at  the  outer  edges  of  the  first  quadrant.  These  conditions 
are  satisfied  by  means  of  a  boundary  collocation  procedure  based  on  the 
displacements  and  the  resultant  forces.  In  the  numerical  computations,  the 
traction- free  conditions  of  the  cracks  are  replaced  by  M  relations,  and 
those  of  the  outer  edges  by  2N  relations. 

The  traction- free  conditions  on  the  cracks  are  common  to  all  considered 
problems.  In  order  to  obtain  M  relations  from  the  traction- free  conditions 
on  crack  A^B^  (stated  in  terms  of  resultant  force  components  P^) ,  we  divide 

this  crack  into  M  unequal  intervals,  D^D2 ,  D2°3 ’ -  DM°M+1 ’  w^ere  °2 ’  D3’*~ 

-,  are  the  mid-points  of  previously  defined  equal  intervals  (Fig. 3). 
Then  the  traction- free  conditions  are: 


Along  crack  A^B^ 


0  (k  -  1,2, - ,M) 


To  impose  the  boundary  conditions  at  the  outer  edges,  EF  and  FG  are 
divided  into  N^  and  N2  equal  intervals,  respectively,  as  shown  in  Fig.  4. 
The  corresponding  interval  lengths  then  are 

hl  "  N  ’  h2  N„  (14) 


The  integer  N  in  Eqn.  (8)  is  chosen  such  that 


N  -  +  N2 


Now  the  boundary  conditions  for  all  five  problems  given  in  Table  1  are 
replaced  by  the  following  relations  in  term  of  the  resultant  force  and  the 


displacement  components : 


Problem  (1) 


s  s 

Along  EF  :  (Px)sk+1-  [Py]sk+1  -  0  (k  -  1,2,---,^) 

k  k 

—  ^k+l  ^k+1 

Along  FG  :  [Pxls  -  0,  [Py]s*+ -  °2h2  <k  -  ^+1,^+2,  — ,N) 
k  k 

Problem  (2) 


(16) 


—  ^k+1  F  , 

Along  EF  :  [u]gK  A-  0,  [PJ*-  ^  c,  [Py]g  -  0  (k  -  1,2,---,^) 

k  Ic 

—  Sk+1  ^k+1 

Along  FG  :  [Px]^  -0,  [Py]g^  -  <?2  h2  (k  -  N^l.^+2, ---,N) 
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Problem  (4) 


Along  EF  :  [u]sk+1-  0,  [ Px ] ^-  ^  c,  [Py]sk+1-  0  (k  -  1,2, --.N^ 
k  lc 

—  ^k+l  ^k+1  G 

Along  FG  :  0,  0,  [P  ]p-^b  (k  -  N^l.^+2 ,  -  -  -  ,N) 


(19) 
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For  convenience,  N2  is  taken  to  be  an  even  integer,  say  2N3 ,  and  the 
following  relations  (20.1)  and  (20.2)  are  used  : 


s  g 

Along  EF  :  [u]  k+1-  0,  [P  ]  k+1-  0  (k  -  1.2.---.N.) 

k  y  Sk  1 


Along  FG  :  [Py]^+1  -  [P [P^1  -  [F,] 

(u)3k+1  -  [u]“-A  .  [v]^1  -  (v] 

k  SN4-k  Sk  \-k 


N3  ”  N2/2,  N4  “2Nl  +  N2+1  (k-Nj+l.Nj+2,  ---  ,^+Nj) 


(20.1) 


For  the  total  forces: 


[Px]E  ‘  [Px]G  "  2°lc  ’  [Py]F  "  ff2b 


(20.2) 


3.  NUMERICAL  ANALYSIS 

3.1.  Dimensionless  stress  intensity  factors 

The  unknown  quantities  K2n>  M2n  (n~  0.1.2, —  ,N-1)  and  (j-1,2,-- 
-  ,M)  are  determined  by  solving  the  simultaneous  relations  consisting  of 
Eqns.  (13)  and  one  of  Eqns.(16)  to  (20)  depending  on  the  considered  problem. 
The  stress  intensity  factor  at  the  crack  tip  is  then  calculated  by 

iri  ** 

Ki  "  fM  T  (21) 

We  consider  the  following  two  cases  of  uniaxial  loading: 


(a)  Compression  in  the  x-direction  :  a -  -cr^,  “  0 

(b)  Tension  in  the  y-direction  :  -  0, 

where  a is  a  positive  reference  stress. 


Calculations  have  been  performed  for  all  possible  combinations  of  the 
following  values  of  the  geometric  parameters:  b/c  -  2,  1.75,  1.5,  1.25,  1, 
1/1.25,  1/1.5,  1/1,75,  1/2;  c/a  -  5,  6,  7,  8,  9,  10;  L/a  -  0.5,  1.,  1.5, 
2.0,  2.5,  3.0,  3.5,  4.0,  4.5,  5.0,  5.5,  6.0,  6.5,  7.0. 

The  numerical  results  are  summarized  in  terms  of  two  dimensionless 
stress  intensity  factors:  one  is  based  on  the  hole  radius;  and  the  other  is 
based  on  the  crack  length.  These  dimensionless  stress  intensity  factors 
are : 
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Here,  superscript  i  denotes  the  problem  number,  and  superscript  x  or  y 
designates  the  loading  direction. 

As  a/c  tends  to  zero,  the  quantities  in  Eqns.  (22)  for  all  five 
problems  should  converge  to  the  results  for  an  infinite  solid  containing  a 
hole  with  edge  cracks.  The  results  for  this  limiting  problem  are  also  given 
in  terms  of  the  following  dimensionless  stress  intensity  factors: 


.(X) 

*»(x) 

I.O 

TJl(x) 

p(x) 

i.o 

i.O  “ 

F.,o  ' 

^(xa)1* 

,(y) 

K(y) 

Ki,o 

p(y)  _ 

K(y) 

i.o 

i.O 

^(jri)1* 

Fa,  0 

KI,0 

ki,o 

i,0  “ 

^(xi)4 

a,0 

aQ(ir  a)* 

(  a.  -  -o_,  a.  -  0  ) 


(  -  0,  a2  -  aQ  ) 


(  ai  ~  -ffn>  ao  “  Xa^  ) 
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where  subscript  x  or  y  designates  the  loading  direction,  and  ^  and  F^  ^ 
denote  the  values  for  biaxial  loadings. 

3.2.  Accuracy  of  numerical  results 

The  accuracy  of  the  numerical  results  would  be  improved  by  increasing 
subdivision  numbers  M,  and  N£.  It  is  however  desirable  to  use  minimum 

values  for  these  quantities,  which  still  provide  needed  accuracy. 

Reasonable  values  of  M  have  been  determined  from  test  calculations  of 
an  infinite  solid.  Table  2  gives  typical  results  for  1/ a -0.4.  As  is 
seen,  any  M  larger  than  3  seems  to  give  reliable  values  with  errors  less 
than  one  per  cent.  Calculations  also  have  been  done  for  several  finite 

regions.  Based  on  this  information,  M  is  taken  to  be  7  to  10  depending  on 
the  region  geometry. 
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Other  series  of  test  calculations  also  have  been  done  to  determine 


Vi 


reasonable  values  for  and  N2.  Tables  3  and  4  give  the  values  of 

p(i,x)  an<j  p(i,y)  f o r  various  values  of  N.  (-N-),  where  b-c,  a/c-0.2, 
cl  1  Z 

i/a -2  and  M-6.  The  convergence  is  excellent,  and  Ng  larger  than  4 
are  likely  to  give  reliable  values  with  errors  less  than  one  per  cent.  On 

the  basis  of  the  above  and  similar  results  for  other  configurations,  and 

N2  are  chosen  to  be  5  to  10  depending  on  the  geometry. 

Thus  the  number  of  unknowns  to  be  determined  from  the  boundary 
conditions  is  confined  to  less  than  50  with  good  accuracy. 


4.  NUMERICAL  RESULTS  AND  DISCUSSION 


4.1.  Infinite  solid  with  a  hole  and  cracks 


Reliable  solutions  are  available  for  an  infinite  solid  containing  a 
hole  with  edge  cracks.  However,  we  have  reanalysed  this  problem  because  of 
its  importance  as  a  common  limiting  case  for  all  other  considered  problems. 
The  results  also  give  useful  information  on  the  accuracy  of  the  basic 
method,  and  help  to  fix  the  value  of  M  to  be  used  for  each  problem. 

The  numerical  values  of  '  and  are  fitted  by  the  following 
five-term  polynomials  for  the  range  0  <  i/a  <  7: 


(a)  Compression  in  the  x-direction 


F^Xq  -  0.0177  -  0.23290  +  1.OO6502  -  1.O75103  +  1.40320*  (0  -  ) 


(24a) 


(b)  Tension  in  the  y-direction 


F^O  “  1-0377  +  0 •  00610  +  2.563302  -  3.O36O03  +  2.78140*  (  0  -  ) 


(24b) 


(c)  Biaxial  loading 


The  results  for  biaxial  loading  are  obtained  by  the  linear  combination 
of  Eqns .  (24a)  and  (24b)  as  follows: 
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(^)"F2,0“ 


(24c) 


Fi.o  -  Fi*0  ■  >fK’  f.,o  -  J  -  Vi  ‘■’I-0’ 


The  numerical  values  from  Eqns.  (24)  based  on  the  present  analysis  are 
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In  excellent  agreement  with  previous  reliable  values  [17]  shown  in  Table  5. 

q  is  usually  convenient  in  theoretical  considerations,  but  in 
analysing  the  test  results,  F  _  would  be  more  useful  since  it  is  normalized 

el  ,  U 

by  the  hole  radius  which  remains  unchanged  during  crack  extension. 

4.2.  Polynomial  expressions  of  stress  intensity  factor 

For  applications,  it  is  desirable  to  represent  the  numerical  results  by 
simple  analytical  expressions.  We  consider  polynomial  expressions  in  terms 
of  the  following  dimensionless  geometric  parameters: 


Our  computations  show  that  the  polynomial  expressions  must  be  such  that 
the  following  conditions  are  satisfied: 

(i)  As  a/c  -*  0,  the  results  for  all  five  problems  converge  to  those 
for  an  infinite  solid,  independent  of  the  values  of  b/c. 

(ii)  Better  accuracy  is  obtained  by  fitting  polynomials  to  F^  than  to 

F  . 
a 

(iii)  For  the  tension  applied  in  the  y-direction  (case  (b)),  the  stress 
intensity  factor  becomes  unbounded  as  a+i-»b. 

The  numerical  values  of  F^i,X^  and  F^i,y^  for  problems  (1)  to  (5)  are 
fitted  by  polynomials  of  the  three  parameters  defined  by  Eqns .  (25).  The 
results  are  summarized  as  follows: 

(a)  Compression  in  the  x-direction  (  -  0  ) 


Case  when  b  >  c 


Case  when  b  <  c 
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(26a. 2) 


(b)  Tension  in  the  y-direction  (  -  0,  ff2  “  a0  ^ 


Case  when  b  >  c 
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(26b. 2) 


(c)  Biaxial  loading 

The  stress  intensity  factor  for  the  biaxial  loading  is  given  by  a 
linear  combination  of  Eqns.  (26a. 1)  to  (26b. 2)  as  follows: 


k£l)  -  (wi)1*  (  -<71F^i’X)  +  a2 F^1,y)  )  (26c) 

Note  that  XpqS  and  YpqS  for  p-0  are  not  included  in  Eqns.  (26a.  1) 
to  (26b. 2)  in  order  to  ensure  that  F^i,x^  and  tend  to  the  results  for 

an  infinite  solid  as  a-a/c->0.  Furthermore,  Eqns.  (26b. 1)  and  (26b. 2) 
are  formulated  in  such  a  manner  that  the  stress  intensity  factor  becomes 
unbounded  as  a  +  £  -*  b . 

Eqns.  (26)  yield  the  dimensionless  stress  intensity  factor  F^ .  To 

obtain  F  from  these  equations,  we  use  the  expression 
£1 


F 

a 


f  (±y 

i'  a  ' 


(27) 


The  coefficients  in  the  polynomials  (26a. 1)  to  (26b. 2)  for  problems  (1) 


to  (5)  are  given  in  Tables  6  to  10,  respectively. 

The  numerical  values  from  these  expressions  are  compared  with  those 
from  direct  analysis  in  Figs.  5  to  9-2,  for  typical  ca s  of  the  five 
treated  problems.  The  agreement  is  quite  good,  and  the  proposed  polynomials 
(26)  are  useful  in  applications. 


5.  SUMMARY  AND  CONCLUSIONS 

A  unified  method  of  analysis  for  various  problems  relating  to  elastic 
plates  containing  circular  holes  with  edge  cracks  is  developed.  The  method 
is  based  on  the  analysis  of  a  unit  rectangular  region  containing  a  circular 
hole  with  edge  cracks,  where  the  boundary  conditions  of  the  outer  edges  are 
suitably  adjusted  in  order  to  treat  various  problems  including  periodic 
arrays  of  holes  with  edge  cracks. 

The  method  is  applied  to  five  important  problems,  and  accurate  values 
of  the  stress  intensity  factors  are  obtained.  These  values  are  fitted  by 
polynomials  which  are  convenient  for  practical  applications. 
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Figure  captions 


Fig.  1-1  Uniform  normal  loads  along  sides  EF  and  FG. 

Fig.  1-2  Strip  with  a  periodic  array  of  holes  and  cracks  in  the 
x-direction. 

Fig.  1-3  Strip  with  a  periodic  array  of  holes  and  cracks  in  the 
y-direction. 

Fig.  1-4  Infinite  solid  with  doubly  periodic  array  of  holes  and  cracks. 
Fig.  1-5  Infinite  solid  with  zig-zag  periodic  array  of  holes  and  cracks 


Fig.  2 


Geometry  of  a  circular  hole  with  edge  cracks. 


Fig.  3 


Subdivision  of  crack  for  numerical  evaluation  of  weighting 


function. 


Fig.  4 
Fig.  5 
Fig.  6 
Fig.  7 
Fig.  8-1 
Fig.  8-2 
Fig.  9-1 
Fig.  9-2 


Subdivision  of  outer  edges  for  numerical  evaluation. 
Ff1,x)  for  Problem  (1),  b/c-2. 

a 

F(2,x)  for  Problem  (2),  b/c-2. 

a 

F(3,y)  for  Problem  (3),  b/c-1. 

F(4,x)  for  Problem  (4),  b/c-1. 

di 

F(4,y)  for  Problem  (4),  b/c-1. 

Si 

F(5,x)  for  Problem  (5),  b/c-2. 
a 


,(5,y) 


for  Problem  (5),  b/c-2. 


M  x} 

Table  3  Variation  of  F'  ’  '  for  the  case  when 

a 

b/c-1 ,  a/c-0.2,  i/a-2,  M-6 


NrN2 

i-1 

i-2 

i-3 

i-4  i-5 

2 

0.2426 

0.0425 

0.0733 

0.0142 

3 

0.2001 

0.0528 

0.0795 

0.0185  j 

4 

0.2156 

0.0514 

0.0793 

0.0182  0.0346 

5 

0.2110 

0.0517 

0.0797 

0.0182  | 

6  ! 

0.2130 

0.0516 

0.0795 

0.0182  i  0.0345 

7 

i  0.2120 

0.0516 

0.0796 

0.0182  ! 

8 

0.2125 

0.0516 

!  0.0795 

0.0182  ;  0.0345 

9  i 

0.2120 

0.0516 

0.0796 

0.0182  , 

10 

0.2119 

0.0517 

0.0795 

0.0182 

P"**  CO  O' 


Table  4  Variation  of  F 


(i,y) 


for  the  case  when 


N1’N2 


i-i 

i-2 

i-3 

i-4 

2.490 

2.191 

_ 

2.265 

2.070 

2.781 

2.239 

2.292 

2.075 

,  2.863 

2.233 

2.294 

2.077 

j  2.829 

2.236 

2.298 

|  2.848 

2.234 

2.295 

2.077 

|  2.837 

2.235 

2.297 

2.077 

2.843 

2.235 

2.296 

2.077 

2.838 

2.235 

2.297 

2.839 

2.235 

2.296 

2.077 

Table  6  Coefficients  of  polynomials  (26.1)  to  (26.2) 


for  Problem  (1)  (  X3qs,  -  0  ) 


p 

q 

s 

yU.l) 

pqs 

xa,2) 

pqs 

y(i.d 

pqs 

y(l,2) 

pqs 

1 

0 

0 

-7.4715E-01 

-2.4215E+00 

5 . 7012E+00 

2 . 2056E+01 

1 

0 

1 

2 . 3737E+00 

3 . 8558E+00 

2 . 6320E+01 

-4 . 3445E+01 

1 

0 

2 

-  5 . 3173E-01 

-1.4332E+00 

-2 . 5258E+01 

1 . 9673E+01 

1 

0 

3 

0 . 0000E+00 

0. 0000E+00 

0.  OOOOE+OO 

0. OOOOE+OO 

1 

1 

0 

-1 . 9973E+00 

1 . 8117E+00 

-1 . 2146E+02 

-  5 . 6748E+01 

1 

1 

1 

-6 . 1593E+00 

-2.0771E+00 

-5.4035E+01 

1 . 2814E+02 

1 

1 

2 

3 . 2767E+00 

6 . 2703E-01 

1 . 3601E+02 

-  5 . 8552E+01 

1 

1 

3 

0.0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

0.  OOOOE+OO 

1 

2 

0 

2.4177E+01 

0. 0000E+00 

3 . 5453E+02 

- 1 . 5351E+02 

1 

2 

1 

-1 . 7526E+01 

0. OOOOE+OO 

5 . 5240E+01 

1 . 9721E+02 

1 

2 

2 

-2 . 9660E+00 

0. OOOOE+OO 

-  3 . 6728E+02 

-8 . 9138E+01 

1 

2 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

1 

3 

0 

-2.0392E+01 

0. OOOOE+OO 

- 1 . 6091E+02 

2 . 7113E+02 

1 

3 

1 

6 . 2492E+00 

0. OOOOE+OO 

-3.0702E+02 

-4 . 1494E+02 

1 

3 

2 

1 . 7132E+01 

0 . OOOOE+OO 

4. 7679E+02 

1 . 9012E+02 

1 

3 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

0 

0 

4 . 5632E+00 

3 . 3112E+01 

5 . 5254E+01 

- 1 . 2963E+02 

2 

0 

1 

-2 . 5129E+00 

-5.0744E+01 

-1 . 8464E+02 

2 . 8193E+02 

2 

0 

2 

-1.4135E+01 

1 . 8632E+01 

2 . 9894E+01 

- 1 . 6561E+02 

2 

0 

3 

0. 0000E+00 

0 . OOOOE+OO 

0. OOOOE+OO 

0 . OOOOE+OO 

2 

1 

0 

-5 . 0766E+00 

-2 . 6862E+01 

4 . 3826E+02 

8 . 1074E+02 

2 

1 

1 

4 . 0657E+01 

4 . 4293E+01 

- 1 . 0866E+02 

- 1 . 5208E+03 

2 

1 

2 

5 . 2748E+01 

-1. 3934E+01 

3 . 7345E+02 

8 . 1247E+02 

2 

1 

3 

0. 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

2 

0 

-1.1439E+02 

0. OOOOE+OO 

-2 . 3193E+03 

-5 . 1824E+02 

2 

2 

1 

8.4083E+01 

0. OOOOE+OO 

2 . 5729E+03 

8 . 2030E+02 

2 

2 

2 

-1 . 2571E+02 

0.  OOOOE+OO 

- 1 . 6524E+03 

-4 . 7128E+02 

2 

2 

3 

0. 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

3 

o 

1 . 4261E+02 

0. OOOOE+OO 

1 . 8022E+03 

-4.2400E+02 

2 

3 

1 

-1.0036E+02 

0. OOOOE+OO 

- 1 . 8598E+03 

9 . 4728E+02 

2 

3 

2 

4 . 1161E+01 

0. OOOOE+OO 

8 . 7650E+02 

-4 . 6230E+02 

2 

3 

3 

0. 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

0 . OOOOE+OO 

Table  7  Coefficients  of  polynomials  (26a. 1)  to  (26b. 2) 
for  Problem  (2) 


p 

q 

s 

y(2 , 1) 

pqs 

v(2 ,2) 
pqs 

y(2,1) 

pqs 

y(2,2) 

pqs 

1 

0 

0 

4.4980E-02 

1 . 7911E-01 

-3.4397E+01 

-9 . 7220E-01 

1 

0 

1 

1 . 4205E+00 

1 . 1863E-02 

1 . 3797E+02 

-6 . 0517E+00 

1 

0 

2 

-1 . 3067E+00 

-4 . 6399E-03 

-1 . 2161E+02 

1 . 2914E+00 

1 

0 

3 

0 . 0000E+00 

0. OOOOE+OO 

1 . 6083E+00 

0. OOOOE+OO 

1 

1 

0 

1.0267E+00 

1 . 0965E+00 

1 . 5796E+02 

1 . 9387E+01 

1 

1 

1 

-8 . 1648E+00 

-  2 . 6270E+00 

-4. 8040E+02 

1.4744E+01 

1 

1 

2 

7 . 6367E+00 

7 . 6423E-01 

4.0362E+02 

-6 . 0552E+00 

1 

1 

3 

0 . 0000E+00 

0. OOOOE+OO 

-2 . 3961E+01 

0.  OOOOE+OO 

1 

2 

0 

-2 . 3842E+01 

-1 . 3997E+00 

-4 . 7568E+02 

-5 . 1005E+01 

1 

2 

1 

6 . 1551E+01 

4 . 7890E+00 

2 . 5405E+02 

-1 . 5115E+01 

1 

2 

2 

-4. 1235E+01 

-1 . 6465E+00 

-2 . 7887E+02 

1 . 1186E+01 

1 

2 

3 

0 . 0000E+00 

0. OOOOE+OO 

5 . 3624E+02 

0. OOOOE+OO 

1 

3 

0 

5 . 0970E+01 

0. OOOOE+OO 

-  2 . 7518E+02 

4 . 2245E+01 

1 

3 

1 

- 1 . 1814E+02 

0. OOOOE+OO 

2 . 7326E+03 

-1 . 1905E+00 

1 

3 

2 

7 . 2459E+01 

0. OOOOE+OO 

- 1 . 3206E+03 

-5 . 5917E+00 

1 

3 

3 

0 . 0000E+00 

0. OOOOE+OO 

- 1 . 3142E+03 

0. OOOOE+OO 

2 

0 

0 

1 . 1656E+00 

-3.4406E+00 

2 . 2052E+02 

0. OOOOE+OO 

2 

0 

1 

-1 . 8971E+01 

2.4326E+00 

-4 . 0462E+02 

0. OOOOE+OO 

2 

0 

2 

1.5030E+01 

-7 . 8577E-01 

-7 . 2028E+02 

0. OOOOE+OO 

2 

0 

3 

0 . OOOOE+OO 

0 . OOOOE+OO 

1 . 0938E+03 

0. OOOOE+OO 

2 

1 

0 

2 . 2888E+01 

1.5467E+01 

-6 . 3733E+02 

0. OOOOE+OO 

2 

1 

1 

6 . 9402E+01 

- 1 . 0844E+01 

1 . 3699E+03 

0. OOOOE+OO 

2 

1 

2 

-8 . 2563E+01 

3.4696E+00 

1 . 2412E+03 

0. OOOOE+OO 

2 

1 

3 

0 . G000E+00 

0 . OOOOE+OO 

-2 . 7180E+03 

0. OOOOE+OO 

2 

2 

0 

-9.0318E-02 

-6.2833E+00 

1 . 8033E+03 

0. OOOOE+OO 

2 

2 

1 

-2 . 4824E+02 

-6 . 1725E+00 

9 . 3307E+02 

0. OOOOE+OO 

2 

2 

2 

2 . 4491E+02 

4.4776E-01 

6 . 8139E+03 

0. OOOOE+OO 

2 

2 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

-9 . 8197E+03 

0. OOOOE+OO 

2 

3 

0 

- 1 . 3510E+02 

0. OOOOE+OO 

4.0440E+03 

0. OOOOE+OO 

2 

3 

1 

4 . 9533E+02 

0. OOOOE+OO 

-1 . 6234E+04 

0 . OOOOE+OO 

2 

3 

2 

3 . 7369E+02 

0. OOOOE+OO 

-2 . 3441E+04 

0. OOOOE+OO 

2 

3 

3 

0. 0000E+00 

0. OOOOE+OO 

3 . 7646E+04 

0 . OOOOE+OO 

3 

0 

o 

0 . OOOOE+OO 

0. OOOOE+OO 

-2.2984E+02 

0. OOOOE+OO 

3 

0 

1 

0. 0000E+00 

0. OOOOE+OO 

-4 . 3422E+02 

0. OOOOE+OO 

3 

0 

2 

0. OOOOE+OO 

0. OOOOE+OO 

2 . 6958E+03 

0. OOOOE+OO 

3 

0 

3 

0. OOOOE+OO 

0. OOOOE+OO 

-2 . 8198E+03 

0. OOOOE+OO 

3 

1 

0 

0 . OOOOE+OO 

0. OOOOE+OO 

-8 . 0843E+02 

0. OOOOE+OO 

3 

1 

1 

0.  OOOOE+OO 

0. OOOOE+OO 

2 . 6784E+03 

0 . OOOOE+OO 

3 

1 

2 

0. OOOOE+OO 

0. OOOOE+OO 

4 . 2135E+03 

0. OOOOE+OO 

3 

1 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

-2 . 3381E+03 

0. OOOOE+OO 

3 

2 

0 

0 . OOOOE+OO 

0. OOOOE+OO 

9 . 7915E+03 

0. OOOOE+OO 

3 

2 

1 

0 . OOOOE+OO 

0. OOOOE+OO 

-4 . 0365E+04 

0. OOOOE+OO 

3 

2 

2 

0. OOOOE+OO 

0 . OOOOE+OO 

-2 . 8534E+04 

0. OOOOE+OO 

3 

2 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

5 . 6549E+04 

0. OOOOE+OO 

3 

3 

0 

0 . OOOOE+OO 

0. OOOOE+OO 

-  2 . 4737E+04 

0. OOOOE+OO 

3 

3 

1 

0. OOOOE+OO 

0. OOOOE+OO 

6 . 3754E+04 

0. OOOOE+OO 

3 

3 

2 

0. OOOOE+OO 

0. OOOOE+OO 

1 . 1662E+05 

0. OOOOE+OO 

3 

3 

3 

0. OOOOE+OO 

0. OOOOE+OO 

- 1 . 5932E+05 

0. OOOOE+OO 

Table  8  Coefficients  of  polynomials  (26a. 1)  to  (26b. 2) 


for  Problem  (3)  (  X__  ,  -  0  ) 

33s  33s 


pqs 

y(3 , i) 
pqs 

y(3 i 2> 
pqs 

y(3,1) 

pqs 

1 

0 

0 

4 . 6011E-03 

8 . 0875E-01 

-1 . 5119E+00 

1 

0 

1 

- 1 . 0581E-01 

-1. 3423E+00 

- 1 . 3759E+01 

1 

0 

2 

1 . 5066E-01 

0.  OOOOE+OO 

9 . 1954E+00 

1 

0 

3 

0 . 0000E+00 

0. OOOOE+OO 

0 . OOOOE+OO 

1 

1 

0 

-  3 . 5730E+00 

-2.0032E+00 

-4 . 9449E+01 

1 

1 

1 

3 . 8074E+00 

4 . 6367E+00 

1 . 5889E+02 

1 

1 

2 

-1 . 2534E-01 

0. OOOOE+OO 

-  8 . 1243E+01 

1 

1 

3 

0. 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

1 

2 

o 

2 . 0181E+01 

5 . 6744E+00 

2 . 4138E+02 

1 

2 

1 

-2 . 9163E+01 

-6 . 9607E+00 

-4 . 7704E+02 

1 

2 

2 

6 . 0356E+00 

0. OOOOE+OO 

1 . 6456E+02 

1 

2 

3 

0. 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

1 

3 

0 

-2 . 5633E+01 

0 . OOOOE+OO 

-2 . 5499E+02 

1 

3 

1 

4 . 5789E+01 

0 . OOOOE+OO 

3 . 6793E+02 

1 

3 

2 

-1.4734E+01 

0 . OOOOE+OO 

-4 . 5176E+01 

1 

3 

3 

0 . 0000E+-00 

0 . OOOOE+OO 

0 . OOOOE+OO 

2 

0 

0 

-3 . 3853E-01 

-1.0764E+01 

1 . 5635E+01 

2 

0 

1 

5 . 3594E+00 

2 . 0515E+01 

5 . 6553E+01 

2 

0 

2 

-  5 . 2904E+00 

0.0000.  K)0 

-3 . 7740E+01 

2 

0 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

1 

0 

8 . 8234E+00 

-1 . 5910E+01 

-1 . 9539E+02 

2 

1 

1 

-4.0061E+01 

-3 . 5513E+01 

2 . 7947E+01 

2 

1 

2 

3 . 7302E+01 

0 . OOOOE+OO 

-  5 . 3576E+01 

2 

1 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

2 

0 

-7 . 1578E+01 

-6 . 9348E+00 

7 . 3457E+02 

2 

2 

1 

1.4890E+02 

4. 9963E+01 

- 1 . 6397E+03 

2 

2 

2 

-8 . 0834E+01 

0 . OOOOE+OO 

1 . 4676E+03 

2 

2 

3 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

3 

0 

9 . 9964E+01 

0. OOOOE+OO 

:  -  9 . 6305E+02 

2 

3 

1 

-2.0570E+02 

0. OOOOE+OO 

S  3 . 0293E+03 

2 

3 

2 

9 . 3005E+01 

0. OOOOE+OO 

;  -2 . 5766E+03 

2 

3 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

;  0. OOOOE+OO 

3 

0 

0 

0. OOOOE+OO 

5 . 2036E+01 

!  0 . OOOOE+OO 

3 

0 

1 

0. OOOOE+OO 

-8 .4383E+01  1  0. OOOOE+OO 

3 

0 

2 

0 . OOOOE+OO 

0. OOOOE+OO 

i  0. OOOOE+OO 

3 

0 

3 

0 . OOOOE+OO 

0. OOOOE+OO 

I  0. OOOOE+OO 

3 

1 

0 

1  0. OOOOE+OO 

5 . 4394E+01 

i  0. OOOOE+OO 

3 

1 

1 

■  0 . OOOOE+OO 

1 . 4518E+02 

;  0. OOOOE+OO 

3 

1 

2 

i  0. OOOOE+OO 

0. OOOOE+OO 

;  0 . OOOOE+OO 

3 

1 

3 

'  0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

3 

2 

0 

j  0 . OOOOE+OO 

i  - 1 . 1444E+02 

0. OOOOE+OO 

3 

2 

1 

0. OOOOE+OO 

1  -  8 . 3128E+01 

0. OOOOE+OO 

3 

2 

2 

0.  OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

3 

2 

3 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

y(3,2) 
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1.0891E+oT 
1.7244E+01 
1.4304E+01 
3 . 9456E+00 
1 . 1683E+01 
8 . 5217E+00 
5 . 4158E+00 
2 . 1868E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
5 . 6818E+01 
1 . 4930E+02 
9 . 7176E+01 
1 . 9697E+01 
5 . 8317E+01 
-2 . 2321E+02 
1 . 6611E+02 
-3.4682E+01 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0.0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . OOOOE+OO 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . OOOOE+OO 
0 . 0000E+00 
0. 0000E+00 


Table  9  Coefficients  of  polynomials  (26a. 1)  to  (26b. 2) 


for  Problem  (4)  (  X^,  X^.  Y^,  Y3qs  -  0  ) 


p 

q 

s 

y(4,D 

pqs 

x(4,2) 

pqs 

y(4,1) 

pqs 

y(4,2) 

pqs 

1 

0 

0 

-  3 . 2587E-01 

1 . 2001E-02 

-6 . 6054E+00 

-4 . 7946E+00 

1 

0 

l 

3 . 7426E-01 

3 . 8127E-02 
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1 

0 

2 
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4 . 8180E-02 
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0. OOOOE+OO 

1 

0 
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0 . 0000E+00 

0. 0000E+00 
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1 

1 

o ! 
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-6 . 7988E-01 
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1 

1 
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-1 . 3825E+00  1 

-  7 . 0045E-01 

1 

1 

2  ! 

1 . 1054E-01 

1 . 6759E-01 

0. OOOOE+OO 

0 . OOOOE+OO 

1 

1 

3  ! 

0 . OOOOE+OO 

0 . 0000E+00 

0. OOOOE+OO 

0 . OOOOE+OO 

1 

2 

0  | 

5 . 8229E-01 

8 . 1387E-01 

-2 . 2990E+01 

-2 . 7635E+01 

1 

2 

1 

1.2307E+00 

2 . 5439E+00 

-4 . 4066E+00 

4 . 2335E+00 

1 

2 

2 

2 . 7677E-01 

-1 . 0634E+00 

0. OOOOE+OO 

0. OOOOE+OO 

1 

2 

3 

0.0000E+00 

0  0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

1 

3 

0 

0 . OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

1 

3 

1 

0 . 0000E+00 

0 . OOOOE+OO 

0. OOOOE+OO 

i  0. OOOOE+OO 

1 

3 

2 

0. 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

1 

3 

3 

0. 0000E+00 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

0 

0 

i  1 . 4991E+00 

-1.4935E+00 

2 . 5284E+01 

1 . 8234E+01 

2 

0 

1 

!  1.0240E+00 

3 . 0513E+00 

-2 . 2889E+00 

1  1 . 6523E+00 

2 

0 

2 

-2 . 7217E+00 

- 1 . 8975E+00 

0. OOOOE+OO 

!  0 . OOOOE+OO 

2 

0 

3 

!  0. 0000E+00 

0. OOOOE+OO 

0 . OOOOE+OO 

0. OOOOE+OO 

2 

1 

0 

- 1 . 2694E+01 

'  6 . 5026E+00 

- 1 . 3826E+02 

■  - 1 . 0653E+02 

2 

1 

1 

9 . 0147E+00 

;  -  6 . 9495E+00 

j  2 . 7114E+01 

;  1 . 4958E+01 

2 

1 

2 

7 . 1809E+00 

j  4 . 6407E+00 

1  0 . OOOOE+OO 

0. OOOOE+OO 

2 

1 

3 

0. 0000E+00 

1  0. OOOOE+OO 

t  0. OOOOE+OO 

i  0. OOOOE+OO 

2 

2 

0 

5 . 1746E+00 

{  -  9 . 5870E+00 

!  1 . 4436E+02 

1  1 . 3627E+02 

2 

2 

1 

-1 . 9962E+00 

9 . 1101E-01 

- 1 . 1742E+01 

'  - 3 . 5588E+01 

2 

2 

2 

- 1 . 3561E+01 

-2 . 9050E+00 

1  0. OOOOE+OO 

:  0. OOOOE+OO 

2 

2 

3 

0. 0000E+00 

1  0. OOOOE+OO 

1  0. OOOOE+OO 

0. OOOOE+OO 
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Table  10  Coefficients  of  polynomials  (26a. 1)  to  (26b. 2) 


for  Problem  (5)  (  X3qs .  Y3qs  -  0  ) 


p 

q 

s 

y  (  5  *  1 ) 

pqs 

X<5’2) 

pqs 

y(5,D 

pqs 

1 

0 

0 

1 . 7531E-01 

3 . 0014E-01 

3 . 0181E+01 

1 

0 

l 

- 1 . 1768E-01 

-6.4930E-01 

-  6 . 3848E+01 

1 

0 

2 

-4 . 6124E-02 

5 . 5971E-01 

2.6805E+01 

1 

0 

3 

0 . 0000E+00 

- 1 . 1408E-01 

0. OOOOE+OO 

1 

1 

0 

-2 . 5995E+00 

2 . 8194E+00 

-2 . 4054E+02 

1 

1 

1 

7 . 0347E+00 

-4 . 2148E+00 

5 . 4821E+02 

1 

1 

2 

-4 . 3277E+00 

8 . 1776E-01 

-  2 . 8360E+02 

1 

1 

3 

0 . OOOOE+OO 

7 . 8482E-02 

0. OOOOE+OO 

1 

2 

o 

0 . OOOOE+OO 

-2 . 9033E+00 

3 . 9446E+02 

1 

2 

1 

0. 0000E+00 

3 . 0314E+00 

-9 . 3436E+02 

1 

2 

2 

0. 0000E+00 

2 . 2003E+00 

5 . 1613E+02 

1 

2 

3 

0. OOOOE+OO 

- 1 . 2209E+00 

0. OOOOE+OO 

1 

3 

0 

0. OOOOE+OO 

0 . OOOOE+OO 

0. OOOOE+OO 

1 

3 

1 

0. OOOOE+OO 

0 . OOOOE+OO 

0 . OOOOE+OO 

1 

3 

2 

0. OOOOE+OO 

0 . OOOOE+OO 

0. OOOOE+OO 

1 

3 

3 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

0 

0 

2 . 5910E+00 

-6 . 2418E-01 

1  -6 . 2190E+01 

2 

0 

1 

-8 . 5914E+00 

-8 . 1785E-01 

!  2 . 1830E+01 

2 

0 

2 

5 . 7373E+00 

1 . 4816E+00 

i  6 . 2948E+01 

2 

0 

3 

0. OOOOE+OO 

-8.6639E-01 

i  0. OOOOE+OO 

2 

1 

0 

8 . 9134E+00 

-  3 . 3216E+00 

1  8.4793E+02 

2 

1 

1 

- 1 . 6126E+01 

8 . 0049E+00 

- 1 . 4123E+03 

2 

1 

2 

6 . 9518E+00 

-  2 . 5632E+00 

4 . 5533E+02 

2 

1 

3 

0.  OOOOE+OO 

1 . 4334E+00 

1  0. OOOOE+OO 

2 

2 

0 

0. OOOOE+OO 

1 . 7599E+01 

!  -1.6084E+03 

2 

2 

1 

0 . OOOOE+OO 

-2.9872E+01 

3 . 1919E+03 

2 

2 

2 

0. OOOOE+OO 

7 . 1522E+00 

'  - 1 . 4621E+03 

2 

2 

3 

0. OOOOE+OO 

I  - 1 . 5332E+00 

0. OOOOE+OO 

2 

3 

0 

0. OOOOE+OO 

I  0. OOOOE+OO 

0. OOOOE+OO 

2 

3 

1 

0. OOOOE+OO 

:  0. OOOOE+OO 

0. OOOOE+OO 

2 

3 

2 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

2 

3 

3 

0. OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 
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-6 . 5174E+00 
2 . 8785E+00 
-2 . 3563E+00 
0 . 0000E+00 
-  7 . 4074E+00 
4 . 1677E+01 
-1 . 4110E+0.1 
0 . 0000E+00 
5 . 5070E+01 
- 1 . 2007E+02 
4 . 5461E+01 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
8 . 3148E+01 
-1 . 1222E+02 
4 . 4637E+01 
0 . 0000E+00 
-3 . 2857E+02 
4 . 1580E+02 
-1 . 6118E+02 
0 . 0000E+00 
1 . 1821E+02 
-5 . 7679E+01 
1 . 7903E+01 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
0 . 0000E+00 
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Figure  1-2 
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CHAPTER  III 


ON  MECHANICS  OF  CRACK  GROWTH  AND  ITS  EFFECTS  ON 
THE  OVERALL  RESPONSE  OF  BRITTLE  POROUS  SOLIDS 

M.  Islda*  and  S.  Nemat-Nasser 

Department  of  Applied  Mechanics  and  Engineering  Sciences 
University  of  California,  San  Diego,  La  Jolla,  CA  92093 

Abstract 

This  paper  is  concerned  with  crack  growth  in  brittle  porous  solids 
under  compression  and  its  effects  on  the  overall  response  of  the  material. 
As  a  mathematical  model,  we  consider  an  elastic  solid  containing  a  zig-zag 
array  of  circular  holes  with  a  pair  of  edge  cracks  (two-dimensional 
problem),  and  solve  this  problem  by  using  a  theory  which  gives  numerical 
results  as  accurate  as  desired.  Based  on  the  analytical  results,  we  discuss 
the  crack  growth  process  and  estimate  the  effective  Young's  moduli  as  well 
as  the  stress-strain  relation  for  porous  solids.  Our  computations  show  that 
the  cracks  emanating  from  the  poles  of  the  circular  holes  extend  in  the 
axial  direction  and  grow  --  in  most  cases  in  a  stable  manner,  but  for 
certain  cases  in  an  unstable  manner  during  an  intermediate  loading  state 
as  the  overall  applied  uniaxial  compression  increases,  reaching  a  certain 

*  Present  address:  Department  of  Mechanical  Engineering,  Kurume  Institute 
of  Technology,  2228  Mukuno  Kamitsumachi ,  Kurume  830,  Japan. 


limiting  maximum  length.  This  maximum  crack  length  strongly  depends  on  the 
ratio  of  the  hole  radius  to  the  hole  spacing  in  the  loading  direction.  The 
effective  Young’s  modulus  in  the  direction  of  the  crack  growth  is  basically 
determined  by  the  initial  porosity,  and  is  little  affected  by  the  crack 
length  or  its  growth  regime,  i.e.,  whether  stable  or  unstable.  We  find  that 
the  overall  axial  stress-strain  curve  remains  monotonic,  exhibiting  no  peak 
stress  or  strain  softening,  as  cracks  extend  in  the  axial  direction  and 
reach  their  limiting  length  with  increasing  axial  stress. 


1 .  INTRODUCTION 


For  a  large  class  of  brittle  materials,  failure  often  involves  the 
formation  and  growth  of  tension  cracks.  Under  farfield  overall  tensile 
forces,  such  tension  cracks  initiate  at  preexisting  flaws  or  defects,  or  at 
inclusions  or  other  material  heterogeneities  which  serve  as  stress 
concentrators.  Tensile  failure  of  this  kind  has  been  extensively  studied 
since  the  pioneering  work  of  Griffith  [1],  It  is  also  known  that  failure  of 
brittle  materials  under  farfield  compressive  forces  is  often  caused  by  the 
formation  and  growth  of  tensile  microcracks  at  preexisting  flaws  or  other 
heterogeneities.  In  this  case,  even  though  the  overall  farfield  stresses  may 
be  compressive,  the  presence  of  defects  may  produce  high  tensile  stresses  at 
the  local  level,  leading  to  the  formation  of  tension  cracks.  Brittle 
failure  of  this  kind  is  rather  complex  and  has  received  attention  of 
researchers  only  recently. 

Various  mathematical  models  have  been  considered  in  order  to  explain 
and  analytically  quantify  brittle  failure  in  compression.  A  noteworthy 
model  is  a  straight  preexisting  crack  inclined  with  respect  to  the  maximum 
overall  farfield  compression,  producing  tension  cracks  at  its  tips  due  to 
the  relative  frictional  sliding  of  its  faces;  see  [2-11].  Recently,  Nemat- 
Nasser  and  Horii  [8,9],  and  Horii  and  Nemat-Nasser  [10,11],  have  analyzed 
the  elasticity  problem  associated  with  this  model  by  considering  a 
preexisting  thin  straight  flaw  (or  a  set  of  such  flaws)  endowed  with  both 
frictional  and  cohesive  resistance,  and  embedded  in  a  linearly  elastic 
solid.  By  examining  the  growth  of  tension  cracks  which  emanate  from  the 
tips  of  such  a  flaw,  as  well  as  from  the  tips  of  an  interacting  set  of  such 
flaws  which  are  suitably  arranged,  these  authors  show  that  the  model  does 
capture  the  essential  features  observed  experimentally,  i.e., 


axial 


splitting  in  uniaxial  compression,  and  faulting  when  the  axial  compression 
is  accompanied  by  some  moderate  confining  pressures.  In  addition,  by 
including  the  possibility  of  the  formation  of  a  plastic  zone  at  the  crack 
tips,  Horii  and  Nemat-Nasser  [9-11]  show  how  the  model  can  be  used  to 
explain  the  process  of  brittle-ductile  transition.  These  authors  also 
report  on  a  set  of  model  experiments  which  seem  to  corroborate  their 
analytical  predictions. 

Microcracks  may  also  nucleate  at  pores  in  a  porous  solid  and  grow  under 
farfield  compressive  stresses.  This  suggests  an  alternative  model  for  the 
analysis  of  failure  in  compression  of  brittle  solids  containing  micro  - 
cavities.  To  be  able  to  examine  the  physical  consequences  of  a  model  of 
this  kind,  one  requires  a  solution  of  the  problem  of  crack  growth  from  the 
walls  of  cavities  under  prescribed  farfield  stresses.  In  the  two- 
dimensional  case,  the  elasticity  problem  of  cracks  extending  from  the  edges 
of  a  circular  hole  has  been  solved  by  Bowie  [12]  for  an  infinitely  extended 
medium.  The  more  general  problem  of  cracks  emanating  from  an  elliptical 
hole  in  a  plane  has  been  treated  by  several  authors  [13-16]  where  some  of 
the  results  have  been  put  in  a  form  convenient  for  design  purposes  [17,18]. 

Recently,  Samrais  and  Ashby  [19]  proposed  to  use  edge-cracks  emanating 
from  circular  (in  two  dimensions)  or  spherical  (in  three  dimensions) 
cavities,  as  a  basic  model  for  the  analysis  of  brittle  failure  in 
compression.  With  the  aid  of  model  experiments,  these  authors  show  that 
cracks  can  grow  from  pores  and  extend  in  the  direction  of  maximum 
compression.  However,  since  an  exact  solution  for  a  finite  body  containing 
a  circular  cavity  with  edge-cracks  does  not  exist,  these  authors  use 
approximate  estimates  based  on  the  beam  theory.  In  particular,  for  a  solid 
containing  periodically  distributed  holes  with  edge-cracks,  Sammis  and  Ashby 


estimate  the  relation  between  the  overall  uniform  axial  stress  and  the  axial 
strain,  arriving  at  a  stress-strain  relation  with  a  peak  stress  followed  by 
a  descending  portion  (i.e.,  "strain  softening").  In  this  model,  the  cracks 
are  assumed  to  initiate  at  the  poles  of  the  circular  holes  and  to  extend 
straight  in  the  direction  of  axial  compression.  Since  approximate 

it 

expressions  based  on  strength-of-material  considerations  are  used,  it  is 
not  clear  whether  an  exact  solution  of  the  associated  elasticity  problem 
would  also  produce  a  peak  stress  followed  by  "strain  softening". 

In  this  paper  we  solve  the  problem  of  a  zig-zag  array  of  circular  holes 
with  a  pair  of  edge  cracks.  In  the  method  used  the  results  can  be  rendered 
as  accurate  as  desired.  Based  on  our  numerical  results,  we  discuss  the 
crack  growth  process,  the  effective  Young’s  moduli,  and  the  stress-strain 
relation  of  porous  solids.  Our  solution  predicts  a  monotonically  increasing 
axial  stress -strain  relation  for  this  model  problem  even  when  the  cracks 
grow  axially  in  an  unstable  manner. 

2.  THEORETICAL  ANALYSIS  OF  ZIG-ZAG  ARRAY  OF  HOLES  WITH  EDGE  CRACKS 

As  a  mathematical  model  of  porous  solids  with  micro-cracks,  we  consider 
an  infinite  solid  with  a  zig-zag  periodic  array  of  holes  with  edge  cracks, 
as  shown  by  Fig.  1.  Due  to  symmetry,  we  solve  this  problem  using  suitably 
symmetric  stress  potentials  defined  on  the  rectangular  region  OEFGO;  see 

*  In  [19],  the  change  in  potential  energy  due  to  cracking  is  calculated  for 
small  crack  lengths  (i.e.,  1  +  L;1  is  used  in  equation  (22)  of  [19]). 

The  results  are  then  used  to  calculate  the  elastic  modulus  of  the  cracked 
body  for  large  crack  lengths.  This  may  be  a  major  cause  of  the  strain 
softening  prediction  in  [19], 


Fig.  2  for  notation  and  definition  of  various  geometric  quantities.  Figs. 
3-1  and  3-2  show  the  two  typical  hole  distributions  corresponding  to 
b/c  -  0.577  and  b/c-1,  which  are  the  equilateral  triangular  array  and  the 
square  array,  respectively. 

We  consider  the  biaxial  loading  case  in  which 


°1  “  * Vo ’  °2  “  ‘Vo  (1) 

where  a ^  and  a ^  are  the  applied  stresses  at  infinity  in  the  x-  and  y- 
directions,  respectively,  and  a ^  is  a  positive  reference  stress. 

The  analytical  procedure  for  the  solution  of  this  problem  is  outlined 
in  Appendix  1.  Here  we  discuss  the  numerical  results. 


The  numerical  values  of  the  stress  intensity  factor  are  given  in 


terms  of  the  following  dimensionless  quantity,  : 

Kt(A. ,A0) 

n  /  \  \  \  *■  •*» 


Let  E  (in  compression)  and  E  (in  tension)  be  the  effective  Young's 
x  y 

moduli  in  the  x-  and  y-directions  for  the  porous  solid  with  cracks.  They 
are  obtained  from  the  displacements,  u  along  EF  and  v  along  FG,  associated 


with  cases  A, 


1,  A„  -  0  and  A.  -  0,  A-  -  -1,  respectively.  Their 


numerical  results  are  given  in  terms  of  the  "modulus-reduction-factors" . 
and  C  ,  which  are  defined  by: 

y  3 

E  E 

cx  “  if-  C  -  f  (3a) 

x  fc0  y  0 

where 

g 

Eq  -  E  (plane  stress),  Eq  -  j (plane  strain)  (3b) 


Note  that  E^  for  compression  is  not  the  same  as  that  for  tension,  once 
cracks  in  the  x-direction  are  developed.  In  this  case,  (compression)  * 


I 

I 


1 


I 

i 

t 


f 

> 


Ey  (uncracked  porous)  >  E^  (tension  with  cracks). 

and  depend  on  Poisson's  ratio.  (Numerical  results  presented  by 
Fig.  8-1  to  Fig.  13  have  been  obtained  for  plane  strain  with  i/  — 0.3.) 


3.  CRACK  GROWTH  PROCESS  IN  BRITTLE  POROUS  SOLIDS  UNDER  COMPRESSION 
3.1.  Uniaxial  Compression  (A^-l,  >2-0) 

The  solid  curves  in  Figs.  4-1  to  4-3  give  the  dimensionless  stress 
intensity  factors  Fj(l,0)  for  various  values  of  b/c  in  the  case  of  uniaxial 
compression.  In  the  range  b/cSl,  which  includes  the  two  typical  cases 
shown  by  Figs.  3-1  and  3-2,  the  value  of  Fj(l,0)  decreases  monotonically 
with  increasing  values  of  L,  i.e.  an  increase  in  the  applied  stress  is 
required  in  order  to  extend  the  crack  (at  constant  critical  stress  intensity 
factor).  Thus  the  crack  extension  regime  is  stable  in  this  range  of  b/c. 
Fj  reduces  to  zero  at  a  certain  value  of  L,  depending  on  the  values  of  b/c 
and  c/a.  This  trend  becomes  more  pronounced  as  b/c  increases. 

For  larger  values  of  b/c,  unstable  crack  growth  may  occur,  as  shown  by 
Fig.  4-3.  We  observe  that,  for  certain  values  of  c/a,  F^  which  is  first  a 
decreasing  function  of  L,  may  begin  to  increase  with  increasing  crack 
length,  L,  after  a  certain  crack  length  is  attained.  However,  since  the 
ascending  portion  of  the  F^ ,  L-curves  are  followed  by  a  descending  portion, 
the  corresponding  unstable  crack  growth  is  expected  to  cease  and  be  followed 
by  a  stable  one. 

Since  is  constant  (  -  K^c  )  during  the  crack  growth,  we  can  define 
the  normalized  applied  stress  by  the  reciprocal  of  F^ ,  that  is 


a 


Klc 


(4) 


a  ,  L-curves  for  b/c-0.577  and  b/c- 1/0.8  are  given  in  Figs.  5-1  and 
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5- 2,  respectively.  In  the  range  b/c  <  1 ,  a  increases  monotonically  with 
increasing  L,  as  shown  in  Fig.  5-1  for  a  typical  case.  Therefore  the  crack 
growth  is  stable  in  this  range  of  b/c.  For  larger  values  of  b/c,  we 
observe  an  intermediate  stage  where  a  decreases  with  increasing  L,  as  shown 
by  Fig.  5-2.  The  crack  growth  behavior  in  this  case  depends  on  the  test 
conditions.  If  the  specimen  is  tested  under  a  displacement-controlled 
machine,  the  a  ,  L-relations  will  be  as  shown  in  Fig.  5-2,  consisting  of  an 
initial  increasing  stage,  an  intermediate  decreasing  stage,  and  a  final 
increasing  stage.  In  load-controlled  tests,  however,  the  crack  will  extend 
from  L  -  to  L  -  Lj  abruptly,  and  then  a  increases  with  increasing  value 
of  L,  following  the  curve  shown  in  Fig.  5-2.  Therefore  in  this  case,  the 
crack  growth  is  temporarily  unstable. 

As  is  seen  in  Figs  5-1  and  5-2,  the  crack  growth,  except  for  the 
unstable  stage  in  Fig.  5-2,  becomes  more  gradual  as  L  increases,  and  L 
approaches  a  certain  limiting  value  Lmax  when  a  tends  to  infinity.  The 
value  of  can  be  determined  graphically  from  diagrams  similar  to  Figs. 

4-1  to  4-3,  by  finding  the  intersection  of  the  F^. -curves  with  the  L-axis. 

The  variation  of  Lmax  with  b/c  and  c/a  is  displayed  in  Figs.  6-1  and 

6- 2,  respectively.  They  show  that  Lmax  varies  almost  linearly  with  both  b/c 
and  c/a,  and,  hence,  can  be  fitted  by  the  linear  relation 

Lnax  -  -  1.025  +0.003  £  +  -  (  -  0.061  +  0.981  £  )  (5a) 

AC  A 

Results  from  Eqn.  (5a)  are  shown  by  dashed  lines  in  Figs.  6-1  and  6-2. 

Eqn.(5a)  may  also  be  expressed  as 

Lnjax  -  -  1.025  +  0.003  -  +  -  (  0.981  -  0.061  -  )  (5b) 

9i  3l  C 
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which  shows  that  Lmax  is  also  linear  in  b/a  for  fixed  values  of  c/a.  In 
Fig.  6-3  we  have  plotted  the  numerical  values  of  Lmax  for  various  c/a, 
taking  b/a  as  the  abscissa.  It  is  seen  that  L^x  is  almost  independent  of 
c/a,  and  can  be  approximated  by  the  following  linear  expression  of  b/a: 

W  -  0-97  -1)  (5c) 

Eqn.  (5c)  is  sufficiently  accurate,  at  least  for  the  range  of  b/a  and  c/a 
shown  in  Figs.  6-1  and  6-2. 

The  above  b/a  -  dependence  of  Ljj,^  in  a  limited  range  of  b/c  can  be 
understood  from  the  behavior  of  the  stress  flow  around  the  crack  tips,  which 
is  mainly  affected  by  the  relative  horizontal  hole  spacing,  b/a,  and  is 
little  disturbed  by  the  relative  vertical  spacing  c/a. 

3.2.  Effect  of  Lateral  Confining  Pressure 

The  lateral  confining  pressure  considerably  reduces  Kj  and,  hence, 

L  .  As  an  example,  F  and  L  for  a  one  percent  lateral  compression  are 
niflx  l  nicix 

given  by  dashed  curves  in  Figs.  4-1  to  4-3.  In  order  to  examine  the  effect 
of  the  lateral  pressure  on  L  ,  calculations  are  performed  for  various 
values  of  X  (  -  ct„/<7,  )  in  the  typical  case  of  b/c  -  0.577,  and  L  is 
evaluated.  The  results  are  plotted  against  X  in  Fig.  7,  together  with  the 
curve  for  an  infinite  solid  (c/a~*®)  obtained  as  a  limiting  case  of  the 
present  problem.  ^max  reduces  to  zero  at  a  certain  value  of  X  depending  on 
c/a.  These  limiting  values  of  X  for  finite  c/a  are  known  as  the  ratio  of 
the  minimum  and  the  maximum  of  the  circumferential  stress  along  the  holes 
without  cracks  which  are  zig-zagly  distributed  [20).  The  calculated  curves 
in  Fig.  7  are  extended  toward  these  limiting  values  by  dashed  lines. 


4.  MODULUS  CHANGE  AND  STRESS -STRAIN  RELATION 


4.1.  Compression  in  the  Crack  Growth  Direction 

Figs.  8-1  to  8-3  show  the  change  in  the  elastic  moduli  due  to  crack 

extension  for  b/c  -  0.577,  1  and  1/0.8,  respectively.  Values  of  the 

modulus  reduction  factor  C^  are  plotted  by  solid  lines  against  the 

normalized  crack  length,  L,  for  various  values  of  the  hole  size  parameter 

c/a.  We  find  that  C^  decreases  for  smaller  values  of  L,  and  tends  to  be 

almost  constant  for  larger  values  of  L.  Note  that  these  curves  are 

terminated  at  points  where  the  corresponding  limiting  maximum  crack  length, 

L  ,  is  attained, 
max 

Since  C^  depends  on  b/c  and  c/a  in  Figs.  8-1  to  8-3,  it  is  of  interest 

to  check  the  relation  between  C  and  the  porosity  f  -  wa2/(2bc) .  In  Fig. 

x  n 

9,  we  have  plotted  C  against  f  for  various  values  of  b/c,  where  L  is  taken 

X  n 

to  be  1.5  for  convenience.  We  find  that  is  not  much  affected  by  b/c. 
Thus  we  can  conclude  that,  for  all  practical  purposes,  the  modulus  reduction 
factor  depends  only  on  f^,  and  is  almost  independent  of  b/c  and  L,  unless 
b/c  is  extremely  small . 

The  limiting  values  of  when  L-^O  are  obtained  by  analysing  a  zig¬ 
zag  array  of  holes  without  cracks  [20] .  Typical  results  are  plotted  in  Fig. 

10,  and  we  find  that  also  strongly  depends  on  f^  and  is  little 

affected  by  b/c.  These  analytical  values  are  fitted  by  the  following 
polynomial  with  good  accuracy: 

<Vl-0  ’  1  •  2  592fH  +  3'050  £H  '  3354fH 


+  -|  -  0.828  f„  +  1.586  fj  +  7.664  f*  ) 

C  H  H  H 

2 

+  -  [  0.310  f  +  0.836  f2  -  11.680  f*  j 

C  H  M  H 


Let's  make  further  comments  on  for  the  two  typical  cases  shown  by 
Figs.  3-1  and  3-2.  As  is  seen  in  Figs.  8-1  to  8-3,  remains  almost 
constant  after  the  crack  length  reaches  a  certain  value  which  seems  to  be 
around  1.5.  Fig.  11  shows  the  values  of  Cx  for  L=0  and  for  L>1.5. 
These  are  the  upper  and  the  lower  bounds  for  C  .  The  values  for  L-0  are 
well  fitted  by  the  following  polynomials: 


^x^/c-o.syy.L-o 
^Cx^b/c— 1 , L-0 


1  -2.977  fu  +  4.910  f2 

n  n 

1  -2.972  f„  +  2.536  f?, 
H  H 


(7a) 

(7b) 


In  the  special  case  of  b/c-1,  ^x^L-0  can  be  a^so  obtained  by 

transformation  of  the  results  for  a  doubly  periodic  array  of  circular  holes 
distributed  with  equal  spacings  in  both  directions  [21]  (see  Appendix  2). 

Now  we  consider  the  stress-strain  relations  for  this  problem.  Let  a 
and  c  be  the  axial  compressive  stress  and  strain,  respectively.  Then  the 
strain  when  the  crack  length  is  L  is  given  by  the  integral 


e 


L 


(8) 


It  is  oberved  from  Figs.  8-1  to  8-3,  that  remains  nearly  constant 
during  the  crack  growth,  in  a  wide  range  of  values  of  b/c.  In  these  cases 
we  have 


x 


(9) 


which  is  written  in  the  dimensionless  form 


s's's's  '  yH  Mrvv<vv  "  M 


where 


V 

V 

V 

y 


aL(7ra)4  i 


°L  "  KIc 


(ID 


L  Kjc  ~  .vic  ‘I,L 

It  is  to  be  emphasized  that  the  linear  relationship  (10)  holds 
approximately  even  for  cases  where  the  crack  growth  is  temporarily  unstable, 
as  is  seen  in  Figs.  5-2  and  8-3  for  b/c  -  1/0.8.  Fig.  12  gives  the 
<7^  ,  relations  for  this  case  of  b/c. 

As  stated  above,  our  results  do  not  show  any  significant  reduction  in 
the  axial  stiffness  of  the  porous  solid  with  axial  crack  growth,  unless  b/c 
is  extremely  small.  Our  numerical  results  produce  monotonically  increasing 
axial  stress  -  strain  relations  without  any  peak  stress,  even  when  axial 
cracks  grow  in  an  unstable  manner  in  the  axial  direction.  The  effect  of 
such  an  axial  crack  growth  on  the  axial  strain  and  axial  stiffness  is  of 
second  order  which  does  not  change  the  character  of  the  corresponding 
stress-strain  relation. 

4.2.  Young*  s  Modulus  Transverse  to  Crack  Growth  Direction 


The  modulus  reduction  factor  C^,  for  tensile  loading,  is  given  by 
dashed  curves  in  Figs.  8-1  to  8-3.  As  is  expected,  decreases 

considerably  with  crack  growth,  and  the  tensile  load-carrying  capacity  of 
the  specimen  reduces  to  zero  when  the  adjacent  cracks  join.  Thus  is  zero 
when  a  +  i  -*  b  or  L -*  (b/c)  (c/a)  - 1 .  The  curves  obtained  by  numerical 
calculation  may  therefore  be  extended  to  their  corresponding  limiting  values 
on  the  L-axis,  as  shown  by  dotted  lines. 

An  abrupt  drop  in  the  value  of  is  observed  in  Figs.  8-1  to  8-3. 

Fig.  13  gives  the  relation  between  C  and  fT,  for  L— 1.5.  C  is  not 
os  y  H  y 
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much  affected  by  the  variation  in  b/c ,  similarly  to  C  .  However,  it  should 


be  noted  that,  unlike  C  , 

x 


C  suffers  a  considerable  reduction  with 

y 


increasing  L.  Therefore,  Fig.  13  should  be  regarded  as  an  example  for 


L— 1.5.  Cy,  for  tensile  loading,  depends  very  strongly  on  L,  as  well  as  on 


fH- 
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APPENDIX  1  :  ANALYSIS  OF  ZIG-ZAG  ARRAY  OF  HOLES  WITH  EDGE  CRACKS 

Analytical  Concept 

Consider  an  infinite  solid  with  a  zig-zag  periodic  array  of  holes  with 
edge  cracks  as  shown  by  Fig.  1.  Due  to  double  symmetry,  we  solve  this 
problem  using  stress  potentials  defined  on  the  unit  rectangular  region  in 
the  first  quadrant,  OEFGO;  see  Fig. 2. 

The  analysis  is  performed  by  superposing  two  stress  states.  The  first 
one  is  represented  by  the  complex  potentials  <£^(z)  and  V>^(z)  having 
singularities  within  the  hole.  The  second  stress  state  is  properly  singular 
at  the  crack  tips,  and  is  realized  by  distributing  suitable  force  doublets 
along  the  cracks . 

The  above  two  stress  states  are  established  in  such  a  manner  that  they 
automatically  satisfy  the  traction- free  conditions  at  the  hole. 

The  boundary  conditions  on  the  crack  and  at  the  outer  edges  of  the  unit 
region  are  given  in  terms  of  stresses  and  displacements.  In  the  present 
analysis,  the  stress  conditions  are  replaced  by  those  of  the  resultant 
forces  in  order  to  increase  the  accuracy  of  the  numerical  results  [22,23]. 

Displacements  and  Resultant  Forces  due  to  the  First  Stress  State 

In  plane  problems  of  elasticity,  components  of  the  stress,  the 
displacement,  and  the  resultant  force  with  respect  to  the  Cartesian 
coordinates  are  expressed  in  terms  of  two  complex  potentials,  4>( z)  and  V>(z). 


as  follows  [24]: 
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2G(u-iv)  -  k^(z)  -  z<£*  (z)  -  V>( z) 


(A1.2) 


Py+  iPx  -  -  <£(z)  -  z^'  (z)  -  0(z) 


(A1.3) 


where  G  is  the  shear  modulus,  and  k  is  defined  in  terms  of  Poisson's  ratio, 


v,  as  follows: 


K  ~  l  +  v  (P^-ane  stress),  -  3  -  kv  (plane  strain)  (Al.4) 

The  complex  potentials  ^(z)  and  ^(z)  for  the  first  stress  state  are 


expressed  as  the  following  Laurent  series : 

*l(z)  -  MQz  +  V[M2n(z2n+1-(2n+l)a4nz'(2n‘1))-K2na4n+2z'(2n+1)] 


^(z)  -  -2MQ  a2z  1  +  Ns:[K  {z2n+1-(2n+l)a4n+4z'(2n+3) 


n-0l  2n 


(Al.5) 


-M  (2n)2a4n+2z‘(2n+1)] 


Using  Eqns. (Al. 5)  in  Eqns . (A1 . 2)  and  (A1.3),  we  have  the  following 


expressions  for  the  displacements  and  the  resultant  forces: 

-  n  n  *  1.  -  .  _  N“X  9n-i.1 


2G(u-iv)(1)  -  ^Q{ 2a2z  +  (/c-l)z}  +  2^  M2n[*z2n+1- (2n+l)zz2n 

+  (2n+l)a4n{-«-(2n’1)-(2n-l)zZ'2n)  +(2n)2  a4n+V  <2n+1>  j 


+  2n  K2n[-2n+1-a4n+2(^'(2n+1)+(2n+l)zz^2n+2>) 
n— 0  *•** 


+(2n+l)a4n+4z'(2n+3)] 


(Al . 6a) 


(Py+1Px)(l)  “  M0<2a2z'1  '2z)+  M2nl-z2n+1-(2n+l)zz2n 

+(2n+l)a4n{z'<2n'1>-(2n-l)zz‘2n)  +<2n)s  a4n+2,- (2n+l) 


+  i  K2n[-z2"n+a4,'+2(2-<2"+1>.<2n*l)2z-(2"+2)| 

n-0  411 


+(2n+l)a4n+V(2n+3)] 


(Al . 6b) 


-  82  - 


The  above  expressions  include  2N  unknown  coefficients  M^n  and  K2 
(n-0,1,2,  —  ,N-1)  to  be  determined  from  the  boundary  conditions. 


Displacements  and  Resultant  Forces  due  to  the  Second  Stress  State 


The  second  stress  state  is  realized  by  distributing  force  doublets  of 
arbitrary  densities  on  the  cracks.  The  force  doublets  are  established  from 
stress  states  for  concentrated  forces  acting  in  an  infinite  solid  which 
contains  a  circular  hole,  and,  therefore,  the  traction-free  conditions  at 
the  hole  edge  are  satisfied  automatically.  Expressions  for  the 
displacements  and  the  resultant  forces  due  to  these  distributed  force 
doublets  are  as  follows  [25] : 

2  i/2 

2G(U-1''><2)  "Si  I 

m-1  -i/2 


s  -  W- 
m  1 


s  -w.. 
m  1 


wl~wl 

(Vwi>’ 


wx  -w2  k(w1-w2) 


-W9  +  (Sni-W,)2  '  (Sm-W2)2  -2(W1-W2)(W1-W2)  z2  (sm-w2)3 


s  -V„  S  -  w0  \ 
m2  m2  m2 


+  2  (  *  z  +—  ) 


z  '  (s  +W_)2 
ra  3 


(Al.7a) 
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where  p  (s  )  (m  — 1,2)  are  unknown  density  functions  of  the  force  doublets, 

mm 

and  the  subscript  m  takes  1  for  the  right  and  2  for  the  left  crack, 
respectively;  the  following  notation  is  also  used  (see  Fig,  Al-  1): 


W,  —  Z  -  X-  ,  W-  —  —  -  X.  ,  W_  -  X- 

1  0,m  2  z  0,m  3  0,m 


X0,m  “  a  +  2  (m"1)  1  x0,m  "  ‘a  *  !  (m_2) 


(Al .8) 


The  density  functions  Pm(sm)  in  Eqns.(Al.7)  are  then  expressed  as  [26] 


P(s  )  -  qm(s  )  [  (|)2  -  S2  ]** 
mm  mm  i  m 


(A1.9) 


It  is  obvious  from  symmetry  that 


q2(-s2)  "  (Al . 10) 

Thus  q^(s^)  is  the  only  unknown  function  in  the  second  stress  state. 

Boundary  Conditions  Based  on  Resultant  Forces  and  Displacements 

A  piecewise  linear  approximation  will  be  used  to  evaluate  the  unknown 
weighting  function  q^(s^).  For  the  right  crack,  this  is  shown  in  Fig.  Al-2, 
where  equally  spaced  intervals  are  used.  Eqns .  (Al.7a)  and  (Al.7b)  now 

reduce  to  linear  expressions  in  f^  (j  -  1,2, — ,M) ,  which  are  the  values  of 
the  weighting  function  at  the  ends  of  these  intervals. 

We  now  have  2N  unknown  coefficients  in  the  first  stress  state  and  M 
unknown  weights,  f j ,  in  the  second  stress  state  to  be  determined  from  the 
boundary  conditions.  Since  both  stress  states  satisfy  the  traction- free 
conditions  on  the  hole,  we  have  only  to  consider  the  boundary  conditions  on 
the  crack  and  at  the  outer  edges  of  the  unit  rectangular  region.  These 


conditions  are  satisfied  by  means  of  a  boundary  collocation  procedure  based 


on  the  displacements  and  the  resultant  forces.  In  the  numerical 
computations,  the  traction-free  conditions  of  the  cracks  are  replaced  by  M 
relations,  and  those  of  the  outer  edges  by  2N  relations. 

Since  the  cracks  are  assumed  to  be  traction- free ,  the  present  analysis 
is  valid  for  such  combinations  of  a ^  and  a ^  that  satisfy  this  assumption. 

In  order  to  obtain  M  relations  from  the  traction- free  conditions  on 
crack  A^B^  (stated  in  terms  of  resultant  force,  P  ) ,  we  divide  this  crack 

into  M  intervals,  D.^ ,  D^D^, —  DmDM+1’  w^ere  D2 ’  °3 ’  —  ’  Dm  are  t^e  m*d‘ 
points  of  the  intervals  shown  in  Fig.  Al-2.  Then  the  traction- free 
conditions  are: 

-  ^k+1 

Along  crack  A  B  :  [P  ]  L  -  0  (k  -  1 ,2 , -  -  - ,M)  (Al.ll) 

y  k 

The  boundary  conditions  at  the  outer  edges  of  the  unit  rectangular 
region  are  given  by  the  following  relations  (A. 12a)  and  (Al.l2b): 


Conditions  for  the  side  EF  : 


u  -  const.  ,  r  -  0 
xy 


(for  any  y) 


Conditions  for  the  side  FG  : 


(a  )  -  (a  ).  ,  (r  )  -  (r  ).  , 

y  x  y  b-x  xy  x  xy  b-x 


-  [v 


tb] 

'b-x 


(Al . 12a) 


(for  any  x) 


Conditions  for  the  total  forces  along  the  outer  edges : 

,G 


[Px]E  '  [Px3G  “  alC 


l*y  3 F  -  ff2b 


(Al . 12b) 


K  H 

regarding  the  relation  [P  ]  -  -  [P  ]  in  Fig.  1. 

X  r  XU 

To  impose  these  boundary  conditions  in  numerical  computations,  EF  and 
FG  are  divided  into  and  N2  equal  intervals,  respectively,  as  shown  in 
Fig.  2,  where  is  taken  to  be  an  even  integer,  say  2N^. 

The  corresponding  lengths  of  the  intervals  on  EF  and  FG  are 


hl  “  2N, 


N„ 


(Al . 13) 


1  ‘'2 

The  boundary  conditions  (A1.12)  are  then  replaced  by  the  following 
relations  (A1.14a)  and  (A1.14b),  in  terms  of  the  resultant  forces  and 
displacements : 

—  ^k+1  ^k+1 

Along  EF  :  [u]s*  -  0,  [P  ]g  -  0  (k  -  1.2,---,^) 

k  ^  k 

~  ,Sk+l  . SN. +l-k  ,sk+l  ,SN,+l-k 

Along  FG  :  |P  1,  -  IP  ]s  4  .  [P  1  -  |P  )  4 

J  k  J  N.  -k  k  N.  -k 

4  4 

r  . Sk+1  .  , SN. +1 -k  .  . Sk+1  .  ,SN.+l-k 
[u]  -  [u]  4  ,  [v]  -  [ v]  4 

k  N.  -k  k  N.  -k 

4  4 

N3  -  N2/2,  N4  -2Nx  +  N2+l  (k-N1+l,N1+2,---,N1+N3) 

(Al . 14a) 
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f  >1 

Iv 

f  'll 


:» 


For  the  total  forces: 

lPx'E  ■  'PJg  -  V  •  lPyl?-"2b  <A1Ub) 

In  order  to  obtain  the  same  number  of  simultaneous  equations  as  there 
are  unknowns,  the  integer  N  in  Eqn.  (Al.6)  is  chosen  such  that 

N  -  +  N2  (A1.15) 

The  unknown  quantities  K2n,  M2n  (n-0,1,2,  —  ,N-1)  and  f^  ( j  —  1,2,  — 

- , M)  are  determined  by  solving  the  simultaneous  equations  (Al.ll)  and 
(A1.14) . 

Accuracy  of  Results 

The  stress  intensity  factor  at  the  crack  tip  is  calculated  by 

Kj  -  fM  Y  4  (A1.16) 

The  numerical  results  are  then  represented  in  terms  of  the  dimensionless 
stress  intensity  factors  given  by  Eqn.  (2). 

The  accuracy  of  the  numerical  results  would  be  improved  with  increasing 
subdivision  numbers  M,  and  N2 .  It  is  however  desirable  to  use  minimum 

values  for  these  quantities  which  still  provide  needed  accuracy. 

Reasonable  values  of  M  have  been  determined  from  test  calculations  of 
an  infinite  solid  containing  a  circular  hole  with  edge  cracks,  as  a  special 
case  of  the  present  problem.  Table  Al-1  gives  typical  results  of  Fj(l,0) 
for  f/a-0.4.  As  is  seen,  any  M  larger  than  5  seems  to  give  reliable 
values  with  errors  less  than  one  percent.  In  computations  of  the  present 
problem,  M  is  taken  to  be  7  to  10  depending  on  the  geometric  parameters. 

Other  series  of  test  c.  lculations  have  been  done  to  determine 
reasonable  values  for  and  N2 .  Table  Al-2  gives  an  example  for  b/c-0.5, 
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c/a- 5  and  i/a- 2.  Fj(l,Q),  C^,  and  are  given  for  various  values  of 
and  N£,  whereas  M  is  fixed  to  be  7.  The  convergence  is  excellent, 
especially  when  -  Nj,  i.e.  h^  -  On  the  basis  of  the  above  and 

similar  results  for  other  configurations,  and  are  chosen  to  be  5  to 
10,  holding  the  relation 

Thus  the  number  of  unknowns  to  be  determined  from  the  boundary 
conditions  is  confined  to  less  than  50  with  good  accuracy. 


1 

I 


APPENDIX  2  :  STIFFNESS  REDUCTION  FACTORS  FOR  ZIG-ZAG  ARRAY  OF  HOLES 

The  computer  program  prepared  in  a  previous  work  [20]  has  been  employed 
to  obtain  the  values  of  (cx)l_0‘  1116  numerical  results  are  then  fitted  by 

Eqn.  (6)  in  terms  of  b/c  and  f H-  Eqns .  (7a)  and  (7b)  are  also  proposed  for 
the  two  typical  arrays  of  holes. 

The  values  for  the  special  case  of  b/c  —  1  can  also  be  obtained  by  the 
transformation  of  the  results  for  the  doubly  periodic  array  of  circular 
holes  distributed  with  equal  spacings  in  both  directions  [21]  (Fig.  A2-1). 
The  procedure  will  be  outlined  below. 

The  stress-strain  relations  for  this  problem  can  be  written  as  follows: 

£U  “  Dllllall  +  D1122a22 

£22  “  D2211<711  +  D2222ff22  (A2.1) 

£12  “  °1212<712 
CT11  "  Cllll£ll  +  C1122£22 
a22  “  C2211£ll  +  C2222£22 
ff12  “  C1212£12 


(A2.2) 


Let  D.  and  C..,in  denote  D..,,  and  C...,  for  the  unperforated  solid 

ijklO  ljklO  ljkl  ljkl 

under  plane  strain.  They  are  given  as  follows: 


1  ^0  1+t/0 

D11110  “  E0  ’  °11220  Eq  ’  D12120  "  EQ 


(A2.5) 


:11110  "  1-*/  2  ’  C11220  “  '  1'0C11110  ’  C12120  ”  1  +  *, 


(A2.6) 


„  E 

O-  1  -  y2  ’  ^0  ’  1  -  „ 


(A2.7) 


where  E,  v  are  Young's  moduluds  and  Poisson's  ratio  of  the  material  under 
plane  stress,  and  E^,  denote  their  equivalent  values  under  plane  strain. 

Values  of  Cllll/ClH10  (  “  fl)  ’  C1122/C11220  (  ”  f2}  and  C1212/C12120 
(  -  f^)  of  this  problem  for  J/-0.3  were  given  by  Nemat-Nasser  et  al.  in  a 

previous  paper  [21],  using  the  notation  C2222/C2222’  C2233/C2233  and  p/p, 

respectively. 

Eqns.  (A2.4)  and  (A2.6)  enable  us  to  calculate  from  f  ’  f2  and  ^3 

by  the  relations 

1  d-^i 
DUU-  E0  q-q 


l  1/Q(  l-vQ2)f2 

h  fl’f2 


(A2.8) 


i  1  +  i/o 


and  the  obtained  values  of  for  ^“0.3  are  given  in  Table  A2-1. 

Now  we  rotate  the  coordinate  axes  by  tt/4,  and  denote  all  quantities  in 
the  rotated  system  of  coordinates  by  superimposed  star.  Then,  the  following 


equations  determine  from 


Table  Al-2. 

Fjd.O).  Cx  and  Cx  for  b/c-0.5,  c/a-5,  i/a-2 ,  M-7, 
where  Neq  (number  of  equations) -M  +  2(N1+N2) 

(plane  strain,  iz-0.3). 


Fj(l,0)  Cx 

0.0343  0.908 
0.0346  0.911 
0.0338  0.891 
0.0345  0.909 
0.0345  0.909 
0.0345  0.908 
0.0345  0.909 
0.0345  0.909 
0.0345  0.908 


0.550 

0.560 

0.484 

0.551 

0.552 

0.548 

0.551 

0.551 

0.551 


3  4  21 

4  4  23 

5  4  25 

5  6  29 

6  6  31 

7  6  33 

7  8  37 

8  8  39 

9  8  41 


1 


Table  A2-1. 


Comparison  of  -  Ex/Eq  for  b/c  -  1 
from  Refs.  [20]  and  [21]  (j/-0.3). 


Cx  for  Zig-zag  holes 


Ref.  [20] 


Fig. 

1 

Fig. 

2 

Fig- 

3-1 

Fig. 

3-2 

Fig. 

4-1 

Fig. 

4-2 

Fig. 

4-3 

Fig. 

5-1 

Fig. 

5-2 

Fig. 

6-1 

Fig. 

6-2 

Fig. 

6-3 

Fig. 

7 

Fig. 

8-1 

Fig- 

8-2 

Fig. 

8-3 

Fig- 

9 

Fig. 

10 

Fig. 

11 

Fig. 

12 

Fig. 

13 

Fig. 

Al-1 

Fig. 

Al-2 

Figure  captions 

Infinite  solid  with  zig-zag  periodic  array  of  holes  and  cracks. 

Unit  rectangular  region  and  notation  of  geometric  quantities. 
Equilateral  triangular  array  of  holes  and  cracks,  b/c  -  0.577. 

Square  array  of  holes  and  cracks,  b/c-1. 

FjCl.O)  and  Fjd.0.01)  for  b/c  -  0 . 577  . 

Fjd.O)  and  Fjd.O.Ol)  for  b/c-1. 

Fj (1 , 0)  and  F^l,  0.01)  for  b/c -1/0.8. 

Relation  between  a  and  L  for  uniaxial  compression,  b/c— 0.577. 

Relation  between  a  and  L  for  uniaxial  compression,  b/c -1/0.8. 

Relation  between  L^x  and  b/c  for  uniaxial  compression. 

Relation  between  Lmav  and  c/a  for  uniaxial  compression. 

Relation  between  L^x  and  b/a  for  uniaxial  compression. 

Relation  between  L^x  and  a2^°\  ^ or  k*axia^  compression,  b/c -0.577 

C  and  C  for  b/c -  0.577. 
x  y  ' 

C  and  C  for  b/c  -  1 . 
x  y 

C  and  C  for  b/c  -  1/0 . 8 . 
x  y 

Relation  between  and  f^  for  zig-zag  holes  with  cracks,  L-1.5. 

Relation  between  C  and  f  for  zig-zag  holes  without  cracks.  (L-0 
x  rl 

Relation  between  C  and  f  for  b/c -  0.577  and  1. 

X  H 

Stress-strain  curves  for  uniaxial  compression,  b/c -1/0.8. 

Relation  between  C  and  f„  for  L-1.5. 

y  H 

Geometry  of  a  circular  hole  with  edge  cracks. 

Subdivision  of  crack  for  numerical  evaluation  of  weighting 


function. 


NOMENCLATURE 


Cartesian  coordinates 
complex  coordinate  -  x  +  iy 
horizontal  and  vertical  spacings  of  holes 
hole  radius 

«a2 

initial  porosity  - 


crack  length 


normalized  crack  length  -  - 

a 


maximum  normalized  crack  length  under  uniaxial  compression 
Young's  modulus  of  material 

Poisson's  ratio  of  material 
3  -  * 


K  — 


^  +  v  (plane  stress) ,  -  3  -  4*  (plane  strain) 
E 


shear  modulus  - 


2(1+*) 


EQ  -  E  (plane  stress)  ,  Eg  -  ^ — ^2  (plane  strain) 


Young's  modulus  of  cracked  porous  solid; 
compression  in  the  x-direction 
Young's  modulus  of  cracked  porous  solid; 
tension  in  the  y-direction 


E 

X 

modulus  reduction  factor  in  the  x-direction  -  — 

E0 

E 

modulus  reduction  factor  in  the  y-direction  - 

E0 


stress  components 

applied  stresses  at  infinity 
positive  reference  stress 


(Aj/O) 


u,  v 


u 

m 

v 

m 

4>(z)  ,  tf(z) 

w 

M 


a 


resultant  force  components 
displacement  components 
mean  value  of  u  along  EF 
mean  value  of  v  along  FG 
complex  stress  potentials 

weighting  functions  of  force  doublet  densities 
subdivision  number  for  the  crack 
number  of  subdivisions  of  EF 
number  of  subdivisions  of  FG 


length  of  intervals  on  EF  - 


2N, 


length  of  intervals  on  FG  -  ^ 


Mode  I  stress  intensity  factor 


K, 


dimensionless  stress  intensity  factor  - 


critical  stress  intensity  factor 


normalized  compressive  stress  - 
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normalized  compressive  strain  - 


-a1  (ira)4  EQ 
KIc  K 
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CHAPTER  IV 

A  MICROCRACK  MODEL  OF  DILATANCY  IN  BRITTLE  MATERIALS 

by 

S.  Nemat-Nasser*  and  M.  Obata** 

ABSTRACT 

For  a  solid  containing  preexisting  flaws,  overall  nonlinear 
constitutive  relations  are  developed  on  the  basis  of  a  model  which  endows  a 
preexisting  flaw  with  frictional  and  cohesive  resistance,  and  which  includes 
nucleation  and  growth  of  tension  cracks  at  the  preexisting  flaw,  as  it 
deforms  under  the  action  of  an  overall  compressive  load.  The  preexisting 
flaws  may  be  randomly  distributed  or  may  have  an  initial  preferential 
distribution.  They  may  be  of  varying  sizes  and  orientations.  Even  when  the 
flaws  are  randomly  distributed,  their  preferential  activation,  and  the 
nucleation  and  growth  of  tension  cracks  at  preferential  flaws,  render  the 
overall  response  of  the  solid  highly  anisotropic.  As  a  first  step  toward  a 
more  complete  constitutive  micromechanical  modeling,  a  dilute  distribution 
of  preexisting  flaws  is  assumed,  rate  constitutive  relations  are  developed 
for  loading  and  unloading,  which  include  hysteresis,  dilatancy,  and  other 
characteristics  observed  experimentally  in  rocks,  ceramics,  concrete,  and 
similar  brittle  materials.  A  number  of  illustrative  examples  are  worked 
out,  and  the  results  are  compared  to  relevant  experimental  observations. 
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1.  INTRODUCTION 


Ceramics,  rocks,  concrete,  and  similar  brittle  materials  have  a  highly 
nonlinear  and  complex  overall  response  to  applied  loads.  This  includes 
load-induced  anisotropy,  hysteresis,  dilatancy,  and  strongly  path- dependent 
stress-strain  relations.  Such  complex  behavior,  to  a  great  extent,  stems 
from  the  activation  of  microdefects  which  are  commonly  present  in  these 
materials.  The  defects  may  be  preexisting  cracks,  cavities,  soft  or  hard 
inclusions,  etc.  These  defects  serve  as  stress  concentrators  and  therefore 
they  locally  change  the  state  of  stress,  leading  to  the  formation  of  tension 
cracks,  even  under  overall  compressive  loads.  It  has  been  demonstrated 
experimentally  and  supported  by  theoretical  models  that  tension  cracks  of 
this  kind  nucleate  and  grow  in  a  highly  preferential  manner  and  therefore, 
even  when  the  microdefects  are  randomly  distributed,  so  that  initially  the 
solid  is  basically  isotropic,  the  response  of  the  solid  becomes  highly 
anisotropic  and  stress-path-dependent  as  loading  continues;  see  Horii  and 
Nemat-Nasser  (1982).  The  failure  modes  in  compression  of  solids  of  this 
kind  have  been  extensively  studied  and  modeled  analytically,  as  well  as 
illustrated  by  model  experiments;  Brace  and  Bombolakis  (1963),  McClintock 
and  Walsh  (1963),  Hoek  and  Bieniawski  (1965),  Scholz  (1968),  Scholz  and 
Kranz  (1974),  Zoback  and  Byerlee  (1975),  Holcomb  (1978),  Holzhausen  (1978), 
Paterson  (1978),  Kachanov  (1982),  Moss  and  Gupta  (1982),  Nemat-Nasser  and 
Horii  (1982),  Kranz  (1983),  and  Horii  and  Nemat-Nasser  (1985a, 1986) .  In 
particular,  it  has  been  shown  by  Nemat-Nasser  and  Horii  (1982)  and  Horii  and 
Nemat-Nasser  (1985a, 1986)  that  axial  splitting,  often  observed  under 
uniaxial  compression,  can  be  explained  in  terms  of  a  model  which  considers 
an  isolated  preexisting  thin  straight  flaw  endowed  with  frictional  and 
cohesive  resistance,  which  nucleates  tension  cracks  at  its  tips  under  the 
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action  of  axial  compression.  The  model  of  a  sliding  crack  with  frictional 
resistance  (but  without  cohesive  resistance)  has  been  criticized  by  some 
authors  (e.g.,  Holcomb  and  Stevens,  1980;  Janach  and  Guex,  1980;  Dey  and 
Vang,  1981;  Costin,  1983)  on  the  grounds  that  actual  electron  microscopic 
observations  do  not  seem  to  support  the  presence  of  such  cracks,  but  rather 
show  a  complex  pattern  of  axial  tension  cracks  emanating  from  a  variety  of 
sources.  There  are,  however,  a  number  of  reasons  which  seem  to  support  the 
usefulness  of  the  basic  model  for  the  micromechanical  constitutive 
formulation  of  the  inelastic  response  and  the  failure  modes  of  materials  of 
this  kind.  These  include;  (1)  The  criticism  of  the  sliding  crack  model  has 
been  based  on  qualitative  aspects  of  microscopic  observations.  (2)  Recent 
data  by  Sondergeld  and  Estey  (1982)  and  Yanagidani  et  al .  (1985),  suggest 
that  the  focal  mechanism  for  stress  induced  acoustic  emission  involves 
double  couples  and  hence  shear  motion.  (3)  Recent  quantitative  stereologic 
evaluation  of  SEM.  analysis  of  rocks  by  Wong  (1985)  shows  that  the  results 
obtained  from  the  sliding  crack  model  are  not  in  conflict  with  SEM 
observations.  Furthermore,  recognizing  the  limitations  of  the  model  (e.g., 
two-dimensionality) ,  Wong  comments  that  the  agreement  between  the 
theoretical  prediction  and  the  quantitative  microscopy  results  is  quite 
remarkable.  (4)  Recently,  Scholz  et  al.  (1986)  have  reexamined  Bridgman's 
ring  experiment  by  subjecting  a  fully  jacketed  ring  of  pyrex  glass  tightly 
fitted  over  a  polished  and  hardened  steel  rod,  to  overall  hydrostatic 
pressure,  and  observed  isolated  axial  cracks  growing  in  a  stable  manner  from 
interior  tiny  flaws  of  no  greater  than  20^m.  Analysis  based  on  Nemat-Nasser 
and  Horii’s  (1982)  sliding  crack  formulation  showed  that  under  the 
prevailing  all-around  compressive  stresses  in  this  experiment,  microflaws  of 
abount  10/im  are  sufficient  to  initiate  such  axial  cracks.  Scholz  et  al. 


(1986)  conclude  that  "axial  cracks  can  be  initiated  by  tiny  flaws  and  grow 
stably  for  long  distances  under  an  overall  compressive  stress  state.  Thus 
shear  cracks  or  other  types  of  stress  concentrations  would  not  necessarily 
be  expected  to  be  prominent  in  SEM  studies.  In  a  heterogeneous  material 
such  as  rock,  where  stress  concentrations  like  those  we  observed  at  the 
ring- rod  interface,  can  be  expected  to  occur  commonly  at  grain  boundaries, 
it  is  not  surprising  that  axial  cracking  predominates."  In  view  of  these 
observations,  and  considering  the  resolution  limit  of  the  SEM  (about 
0.02^m),  criticism  of  the  sliding  crack  model  based  on  the  qualitative 
examination  of  SEM  results  seems  to  require  reevaluation. 

Notwithstanding  the  above  comments ,  we  do  not  intend  to  suggest  the 
sliding  crack  mechanism  as  the  only  source,  or  even  the  major  source,  of 
axial  cracking.  (Indeed,  model  studies  (Nemat-Nasser ,  1985)  have  shown 
raicrocracking  under  axial  compression  from  inclusions  and  cavities  of 
various  geometric  shapes.)  What  we  intend  to  do  is  to  consider  the  sliding 
crack  model  as  a  representative  of  a  thin  flaw  which  can  deform  plastically 
(and  hence  has  cohesive  resistance)  and  can  slide  (and  hence  has  frictional 
resistance),  and  in  a  systematic  and  rigorous  manner  examine  the  results. 
These  results  then  can  be  used  to  study  the  merits  or  shortcomings  of  the 
basic  model. 

To  this  end,  we  shall  consider  a  solid  which  contains  a  dilute 
distribution  of  such  ideal  flaws.  The  flaws  may  be  of  varying  sizes  and 
orientations  which  may  have  either  random  or  preferential  distributions. 
The  matrix  is  assumed  to  be  isotropically  elastic.  Inelasticity  and 
anisotropy  develop  because  of  the  frictional  and  cohesive  resistance  of  the 
preexisting  flaws,  and  because  of  the  formation  and  growth  of  tension  cracks 


at  the  tips  of  the  flaws.  During  unloading,  the  effects  of  frictional 
locking  and  the  cohesive  resistance  of  the  flaws,  their  backsliding  upon 
further  unloading,  the  partial  closure  of  tension  cracks,  and  finally  the 
residual  strains  are  included.  All  calculations  are  made  for  a  two- 
dimensional  model.  Once  the  basic  parameters  of  the  model  are  defined,  then 
the  problem  is  analyzed  consistently  and  systematically,  without  introducing 
additional  assumptions.  This  allows  us  to  make  a  judgment  on  the  merits  of 
the  model  from  the  obtained  results  which  are  then  a  true  reflection  of  the 
basic  model's  assumptions. 


2.  FORMULATION 

We  consider  a  micromechanical  modeling  of  the  inelastic  response  of 
brittle  materials  such  as  ceramics  and  rocks  at  relatively  low  temperatures 
and  under  moderate  confining  pressures  so  that:  (1)  rate  effects  can  be 
ignored;  and  (2)  micro -cracking  and  inelastic  slip  at  preexisting  flaws  can 
be  regarded  as  the  basic  micromechanisms  giving  rise  to  inelasticity;  for 
rather  detailed  discussions  and  extensive  lists  of  relevant  references,  see 
Paterson  (1978),  Kranz  (1983),  and  Horii  and  Nemat-Nasser  (1986).  The  basic 
building  block  of  the  model  is  the  mechanism  shown  in  Fig.  1.  It  consists 
of  an  isolated  preexisting  flaw  PP'  and  tension  cracks  PQ  and  P'Q'  embedded 
in  an  isotropically  elastic  unbounded  solid.  Under  the  action  of  farfield 
compressive  stresses,  the  preexisting  flaw  PP'  undergoes  frictional  sliding 
which  may  be  accompanied  by  cohesive  (due  to  possible  plastic  deformation  of 
the  flaw)  resistance  as  well  as,  possibly,  dilatancy  normal  to  PP' ,  due  to 
the  deformation  of  the  flaw  and  the  possible  presence  of  asperities.  The 
relative  deformation  of  the  two  faces  of  the  flaw  results  in  the  creation  of 
high  tensile  fields  at  the  flaw's  tips  and  leads  to  the  formation  of  tension 
cracks  which  grow  toward  the  direction  of  maximum  compression.  The  model  of 
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frictional  sliding  of  preexisting  cracks  (without  cohesion  or  dilatancy) 
which  leads  to  the  formation  of  tension  cracks  at  its  tips  was  originally 
suggested  by  Brace  and  Bombolakis  (1963)  and  has  been  studied  both  by  model 
experiments  and  analytically  by  a  number  of  investigators;  McClintock  and 
Walsh  (1963),  Holzhauzen  (1978),  Moss  and  Gupta  (1982),  Kachanov  (1982), 
Nemat-Nasser  and  Horii  (1982),  Horii  and  Nemat-Nasser  (1985a ,b , 1986) ,  and 
Steif  (1984) .  Although  micromechanical  observations  on  actual  rocks  have 
not  clearly  identified  preexisting  cracks  as  the  basic  source  of 
inelasticity  in  rocks,  recent  analysis  by  Wong  (1985)  seems  to  suggest  that 
the  model  is  indeed  viable  and  does  capture  both  qualitatively  and 
quantitatively  some  of  the  main  features  of  the  problem. 

In  this  section  we  shall  calculate  the  inelastic  strain  due  to:  (1) 
slip  along  the  preexisting  flaw,  (2)  dilatancy  normal  to  the  flaw,  and  (3) 
the  associated  out-of-plane  tension  cracks  at  the  tips  of  the  flaw.  ^T&l 
dilatancy  normal  to  the  flaw  will  be  assumed  to  be  very  small,  but  its 
effect  on  creating  tension  cracks  at  the  tips  of  the  flaw  is  significant. 

The  total  strain  at  each  state  depends  on  the  parameters  that  define 
the  basic  model  shown  in  Fig.  1.  The  values  of  these  parameters  are  highly 
stress -path-dependent.  Hence,  an  incremental  solution  is  necessary.  Using 
the  computed  strain  tensor,  we  formulate  a  rate -constitutive  relation,  using 
the  equilibrium  conditions  for  a  representative  flaw  and  the  associated 
tension  cracks. 

Throughout  this  paper  we  shall  use  a  fixed  rectangular  Cartesian 
coordinate  system,  xlt  x2,  as  the  background  frame  of  reference.  To 
simplify  notation,  we  use  the  following  direct  representation: 


a:b  “  aijbij> 

“  ^ijklhkl' 

a®b  -  a^b^,  i,j,k,l  -  1,2, 
where  repeated  Indices  are  summed. 

2.1  Description  of  Strain 


The  total  strain 

consists 

of  an  elastic  strain  due 

to 

the 

elastic 

deformation  of  the 

matrix, 

plus  an  inelastic  strain 

due 

to 

slip  and 

dilatancy  of  the  preexisting  flaws,  accompanied  by  inelastic  strain  due  to 
micro -cracking.  In  the  sequel  all  physical  quantities  with  the  dimension  of 
length  are  normalized  by  dividing  by  half  of  the  preexisting  flaw  length,  c. 
For  example,  if  b  is  the  slip  (assumed  constant)  along  flaw  PP' ,  we  use 
b  -  b/c  to  represent  the  dimensionless  slip.  Similarly,  if  d  is  the 
displacement  (assumed  constant  over  the  flaw  length)  of  a  flaw  normal  to  its 
face,  we  use  d  -  d/c. 

Let 

f  -  c2N  (2.2) 

represent  the  density  of  the  flaws  of  length  2c  and  orientation  <j> ,  where  N 
is  the  number  of  such  flaws  per  unit  area.  The  strain  due  only  to  slip  b  of 
these  flaws  then  is 


£b  "  2fbEo- 


where 


-sin  2^  cos  2 <f> 
cos  2 <(>  sin  2$ 


(2.3a) 


(2.3b) 
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Similarly,  the  strain  due  to  dilatancy  d  normal  to  PP'  is  given  by 


£d  "  2fdEl 


(2.4a) 


where 


I 


$ 

$ 


cos2<p  cos  ^sin  ^ 
Hi  cos  <£sin  ^  sin2^ 


(2.4b) 


see  Appendix  A.  Since  the  sliding- induced  opening  of  a  preexisting  flaw 
normal  to  the  surface  must  be  bounded,  we  shall  assume  that  this  opening 
becomes  saturated  and  hence  use 


$ 

‘W 

4 


d  -  d0(l  -  e  iU). 


(2.5) 


Note  that  b  (and  hence  d)  is  assumed  to  be  uniform  along  the  flaw.  This 
simplification  is  reasonable  once  tension  cracks  are  suitably  long. 
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Slip  b  and  dilatancy  d  introduce  normal  and  tangential  concentrated 
gaps  of  values  (b  sin  9  +  d  cos  9)  and  (b  cos  9  -  d  sin  9),  at  points  P  and 
P’ ,  for  cracks  PQ  and  P'Q' ,  of  common  length  i,  and  orientation  9  with 
respect  to  the  preexisting  flaw  PP' .  The  strain  associated  with  these  gaps 
can  be  computed  and  is  given  by 


“  fbiq0  +  fdigj , 


where 


1  -2sin  ^cos (9  +  ^)  cos (9  +  2 4>) 

2  cos(0  +  2<£)  2cos  $sin(0  +  4>) 


St  "T 


2cos  <f>cos(9  +  i i>)  sin(0  +  2 <j>) 

sin(0  +  2 <f>)  2sin  (£sin(0  +  4>) 


(2.6a) 


(2.6b) 


(2.6c) 


see  Appendix  B. 


The  applied  stresses  change  the  crack  opening  displacement.  This  is 
not  included  in  expression  (2.6a).  Under  the  condition  of  a  positive  Mode  I 
stress  intensity  factor  at  Q  and  Q'  (crack  closure  will  be  considered  later 


on),  i.e.,  for  Kx  >  0,  the  strain  due  to  the  change  in  crack  opening  caused 
by  the  applied  loads  can  be  computed  and  is  given  as  follows: 


£a  “  *l2[Mg:a)a  +  <*:£)£], 


(2.7a) 


where 


cos  2(6  +  <f>)  cos(0  +  0)sin(0  +  <j>) 

sin(0  +  <f>)cos(6  +  <f>)  sin2(0  +  4>) 


(2.7b) 


-sin2(0  +  4>)  cos2(0  +  <j>) 
cos2  (8  +  <(>)  sin2  (6  +  <t> )  ’ 


(2.7c) 


see  Appendix  C.  In  (2.7a),  /i  is  the  shear  modulus,  u  is  Poisson's  ratio, 


and  a  is  the  stress  tensor. 


With  L  denoting  the  elastic  modulus  tensor  of  the  matrix  material,  the 
total  strain  tensor  is  now  given  by 

t  —  L' 1 : a 

+  2fbp0  +  2fdpx 

+  fbi!q0  +  fdiqj  (2’8) 

+  f  1  mi2[4(o:g)g  +  (o  :£)§]. 


2.2  Rate  Constitutive  Relations 

To  complete  the  formulation,  we  must  relate  the  kinematic  parameters  b, 
d,  2  and  9  to  the  applied  stress  a.  We  shall  do  this  by  calculating  the 
opening  mode  stress  intensity  factor,  Kt,  at  crack  tips  Q  and  Q'  in  two 


different  ways:  (1)  by  considering  the  applied  loads  and  the  condition  of 
slip  across  the  flaw  (see  Fig.  2a),  as  discussed  by  Horii  and  Nemat- 
Nasser  (1936),  and  (2)  by  calculating  the  stress  intensity  factor  in  the 
presence  of  the  applied  stress  and  the  imposed  concentrated  gaps  at  P  and 
P'  (see  Fig.  2b). 


As  has  been  shown  and  verified  by  numerical  examples  by  Horii  and 
Nemat-Nasser  (1986) ,  the  Mode  I  stress  intensity  factor  at  Q  and  Q'  of 
cracks  PQ  and  P'Q'  can  be  computed  with  very  good  accuracy  by  considering  an 
equivalent  crack  of  length  2i,  subjected  to  a  pair  of  colinear  concentrated 
forces,  F,  as  well  as  the  applied  overall  stresses,  as  shown  in  Fig.  2.  In 
this  representation,  F  denotes  the  driving  force  transmitted  across  the 
preexisting  flaw;  it  lies  along  the  flaw  and  equals  in  magnitude  the 
resultant  force  transmitted  across  the  flaw.  Following  the  procedure 
outlined  by  Horii  and  Nemat-Nasser  (1986),  we  obtain 


2 

y*(i  +  £*) 


(o:g2 


rc)  sin  6  +  a  :a  Jn£  , 


(2.9a) 


where 


P2“Po+MEi.  **-0.27,  rc  “  e'^b^  ,  |b|  *  J^|b|dt.  (2.9b-d) 

Here,  /I  is  the  coefficient  of  friction,  rc  is  the  value  of  cohesion  at  the 
current  value  of  the  total  absolute  value  of  slip  (we  shall  assume  that 
cohesion  is  an  exponentially  decaying  function  of  the  accumulated  slip  with 
decaying  factor  rj  defined  by  (2.9d)),  and  is  the  initial  cohesion.  The 
parameter  t*  is  introduced  in  such  a  manner  that  (2.9a)  is  valid  even  at 
crack  initiation,  i.e.,  at  I  -  0+.  Since  the  opening  mode  stress  intensity 


factor  for  a  closed  sliding  crack  is  maximum  in  a  di.ection  which  makes  an 
angle  9C  -  0.39*  with  the  crack  orientation,  l*  is  obtained  by  equating  the 


stress  intensity  factor  given  by  (2.9a)  at  2  —  0  and  8  —  6C  to  that  maximum 
value.  In  the  present  application  the  effect  of  tension  cracks  is  important 
once  they  become  relatively  long.  Hence,  whether  or  not  2*  is  included  in 
Eq.  (2.9a)  does  not  affect  the  results,  although  its  presence  renders  Kt 
non- singular  at  i  -  0 . 

The  Mode  I  stress  intensity  factor  associated  with  the  gaps  and  in  the 
presence  of  applied  stresses,  on  the  other  hand,  is  given  by 

Kj'  -  ~  <b  sin  6  +  d  cos  * )  +  ’  (2.10a) 

where 

r*  --£**■,  (2.iob) 

see  Appendix  D. 

To  obtain  a  relation  between  b  and  a,  we  require 

Kj  -  Kj '  when  CT:g2  -  tc  >  0.  (2.11) 

Furthermore  the  crack  length  i  is  obtained  such  that 

Kr  -  Kc,  (2.12) 

where  Kc  is  the  critical  value  of  the  stress  intensity  factor  in  the  opening 
mode.  9  is  determined  so  as  to  maximize  Kx  for  a  given  length  2. 

It  should  be  noted  that  model  experiments  (Nemat-Nasser  and  Horii, 
1982)  and  exact  analysis  (Horii  and  Nemat-Nasser,  1983,  1985a)  show  that 

O 

tension  cracks  emanate  from  the  flaw  tips  at  an  angle  of  about  70  with 
respect  to  the  orientation  of  the  flaw,  and  curve  toward  the  direction  of 
the  maximum  compression.  The  model  shown  in  Fig.  1  substitutes  the  curved 


crack  by  an  equivalent  straight  one;  the  equivalence  is  in  the  sense  that 
the  orientation  6  of  the  equivalent  straight  crack  is  computed  at  each 
incremental  loading  by  maximizing  Kx  with  respect  to  6. 

Because  of  the  complex  and  nonlinear  relations  that  exist  among  various 
parameters  which  defined  the  basic  model  shown  in  Fig.  1,  the  overall 
stress-strain  relation  will  be  highly  nonlinear  and  stress-path-dependent. 
It  is  therefore  necessary  to  formulate  the  problem  incrementally.  To  this 
end  we  will  express  l ,  9,  and  b,  in  terms  of  a,  for  given  loading  regimes. 
These  quantities  can  be  written  symbolically  as 

i  -  A :a,  d  -  B:<7,  b  -  C:d.  (2.13a-c) 

We  will  calculate  the  second- order  tensors  A,  B,  and  C  for  prescribed 
loading  regimes.  Note  that,  in  view  of  assumption  (2.5),  we  have 

d  -  f(d0  -  d)b.  (2.14) 

With  Eqs.  (2.13)  and  (2.14),  and  upon  time  differentiation  of  (2.8),  we 
obtain  the  following  rate  constitutive  equation: 

e  -  D:d , 

D  -  L'1 

+  f{bq0  +  dq!  +  — —  -  7ri[4(o:a)a  +  (£:^)^])®A 
+  f{bir0  +  dirj  +  1  8~  ~  *22[U(o:p)a  +  4(a:a)^  +  (o:y)P  + 

+  f  { 2g0  +  iq0  +  C(d0  ’  d)  (2gj  +  iq^jSC 
+  f -^-r — —n22[  4a  ®  a  + 

o  LI  —  —  —  — 


where 


2.3  Specific  Cases  During  Loading 


In  what  follows,  we  consider  various  loading  cases  separately.  We 
shall  use  a  symbolic  representation,  and  give  the  details  of  expressions 
such  as  3Kj/3b  and  bx  in  Appendices  E  and  F. 


(1)  Sliding  occurs  without  crack  growth :  This  happens  when  the  driving 
shear  stress  on  the  flaw  exceeds  the  frictional  and  cohesive  resistance  of 
the  flaw,  i.e.,  when  £:g2  -  rc  >  0,  while  Kx  <  Kc.  From  condition  Kr  -  Kx' , 
we  obtain,  by  time  differentiation, 

r^Kj  aKj'-l  . 

[3b  3b  J  b  “  (2.16) 

where  bx  and  b2  are  the  matrix  coefficients  of  o,  obtained  from  (2.9a)  and 
from  (2.10a),  respectively.  Thus  we  have 


b 


and,  hence, 


£Kj  3K:' 
~3b~  *  ItT 


^2) :£> 


(2.17) 


£  ~  ~JiT  1  gKi>  (b:  -  b2),  A  -  0,  B  -  0.  (2 . 18a-c) 
lb*  ‘  "3b" 

(ii)  Sliding  occurs  with  crack  growth :  This  happens  when  Kr  -  K.c  and 
£:E2  *  rc  -  The  problem  can  be  stated  as  follows:  Maximize  Kx  with 
respect  to  9  and  b  under  the  condition  that  Kx  -  Kf.  With  the  aid  of  a 
Lagrangian  multiplier,  we  arrive  at  the  following  conditions: 


Kx  -  K{  -  Kc, 


1 


'y* 

•:;.i 

;vt 


aKjaKj'  aKjaKj' 

39  3b  ’  ab  39 


j  -  o. 


Time  differentiation  yields , 


3KT  .  dKj  .  3KT  . 

TT'  'St*  +lb  b  -  S.-2. 
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aKt'.  aKt'.  3kt 

TTi+17'  +ltTb-^. 


3J  ,  .  ;  .  j9J_h  h  • 

Ui  +  39  6  +lbh  ~  ~3--' 


(2 . 19a, b) 


(2 . 20a-c) 


These  linear  equations  yield  i,  9,  and  b.  The  solution  is  in  the  form 
(2.13). 


2.4  Specific  Cases  During  Unloading 

After  a  continued  loading  regime,  we  consider  unloading  by  reducing  the 
applied  stresses.  Unloading  at  specific  flaws  may  also  take  place  when,  for 
example,  we  change  the  load  path.  Our  analysis  of  unloading  includes  all 
cases  when  the  net  shear  traction,  ct:£2  -  rc,  transmitted  across  a 
preexisting  flaw  decreases  for  whatever  reason. 

(i)  Crack  grows  without  sliding:  Upon  unloading,  cracks  may  continue  to 
grow  as  the  lateral  compression  is  reduced  but,  because  of  an  initial 
locking  of  the  flaw,  the  maximum  gaps  at  P  and  P'  remain  fixed.  The 
conditions  for  crack  growth  then  are 

3K,' 

Kx'  -  Kc,  -jf  -  0.  (2 . 21a, b) 

Equation  (2.21b)  ensures  that  cracks  grow  in  a  direction  which  maximizes 
Kj'  .  Time  differentiation  now  yields 


v-.t'v 


3Kt'  . 
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(2 . 23a-b) 


From  these  we  have  £  and  6  in  the  form  (2.13).  Since  b  remains  fixed,  in 


this  case  b  —  0 ,  and  hence  C  —  0 . 


(ii)  Backsliding  occurs  without  crack  closure :  When  the  residual  shear 


stress  on  the  preexisting  flaw  overcomes  the  cohesive  and  frictional 


resistance,  backsliding  of  the  flaw  may  take  place.  Since  the  frictional 


and  cohesive  resistance  on  the  preexisting  flaw  acts  against  the  direction 


of  the  relative  displacement,  Kr  in  (2.9a)  must  be  redefined  for  the 


unloading  case  by 


Kl  ”  Jn(£  +  O  +  'c)sin  8  +  o:aJn£, 


(2.24a) 


E2  “  Eo  •  MBi. 


(2.24b) 


Before  backsliding  takes  place,  Kj  of  (2.24a)  may  formally  be  larger  than 


the  actual  stress  intensity  factor  obtained  on  the  basis  of  the  applied 


overall  loads  and  the  existing  gaps  at  P  and  P' .  This  is  similar  to  the 


loading  case  with  o:g2  *  rc  -  0,  so  that  no  slip  occurs.  The  quantity  Kt, 


defined  by  (2.24a),  attains  its  usual  physical  meaning  only  when  backsliding 


occurs,  or  is  just  about  to  occur.  Therefore,  the  condition  for  the 


initiation  of  backsliding  is  given  by 


Kj  <  Kf; 


(2.25) 


see  Appendix  G.  During  backsliding,  however,  we  require  that 


Upon  time  differentiation  at  fixed  i  and  9  we  obtain  an  equation  similar  to 
(2.17)  for  b. 


(iii)  Crack  closure  without  backsliding:  When  crack  closure  occurs,  the 
strain  tensor  is  given  by 


c  —  L‘ 1  :<j 

+  2fbg0  +  2fdg1 


-  f  Tp — z — ^ 7  - ~r  (b  sin  9  +  d  cos  9)2 a 

2jt  1  -  v  (o:a)  - 

im 

+  f  (b  cos  9  -  d  sin  0)/ 3 

+  fls/ 


(2.27) 


where  denotes  the  maximum  crack  length;  see  Appendix  H.  From  the  time 
differentiation  of  (2.27)  and  Eqs.  (2.13),  we  obtain  the  rate  consti£utiv^ 
equation  for  this  special  case  in  the  following  explicit  form: 


c  —  D:<7, 

D  -  L' 1 

+  f  [2g0  +  2f(d0  -  d)gx 

'77^7  ®  +  ’  d)c°s  9)( b  sin  9  +  d  cos^)a8^ 

+  ~2~  {cos  9  -  f(d0  -  d)sin  0  }0]  ®  C 

+  f  l~2^T~^~V  (ct  :a)"2  (b  sin  9  +  d  cos  n£mie>i^ 

Note  that  0  is  fixed  in  this  case. 


The  condition  for  crack  closure  without  backsliding  of  the  flaw  is 


-  138  - 


whose  time  differentiation  at  constant  8  and  b  yields 


aKt' 

*  -  (b2:a)/  (jy  )  -  k:a,  B  -  C  -  0. 


(2. 30a-c) 


(iv)  Crack  closure  with  backsliding:  The  conditions  are 


Kr  -  Kj '  -  0.  (2 . 32a, b) 

Since  8  remains  constant  during  crack  closure,  only  i  and  b  must  be 
expressed  in  terms  of  a.  This  is  done  by  time  differentiation  of  Kj  -  0  and 
Kr'  —  0  at  constant  8  and  by  solving  the  resulting  linear  equations  for  2 
and  b.  B  -  0  because  8  remains  constant. 

2.4  Comments  on  the  Model 

The  model  presented  above  includes  a  number  of  features  which  have  been 
suggested  by  several  researchers  in  the  past,  e.g.  Brace  and 
Bombolakis  (1963),  Zoback  and  Byerlee  (1975),  Kachanov  (1982),  and  Moss  and 
Gupta  (1982).  In  particular,  Moss  and  Gupta  have  obtained  nonlinear 
stress -strain  and  dilatancy  relations  based  on  a  sliding  preexisting  crack 
with  tension  cracks  (wing  cracks)  emanating  from  its  tips.  These  authors 
include  frictional  and  cohesive  effects,  and  calculate  the  inelastic  strain 
associated  with  the  slip  of  the  preexisting  crack  as  well  as  of  the  tension 
cracks.  In  their  model,  however,  the  orientation  of  the  tension  cracks 
(wings)  relative  to  the  preexisting  crack  is  fixed  throughout  the  entire 
loading  history  and  must  be  prescribed.  Moreover,  in  their  model  it  is 
necessary  to  prescribe  an  initial  length  and  crack  opening  displacement  for 
the  wings.  In  the  model  presented  here,  the  flaws  at  which  tension  cracks 
nucleate,  the  orientation  and  the  length  of  the  tension  cracks,  their 


opening  displacement,  and  whether  they  continue  to  grow,  become  dormant,  or 
actually  close,  are  all  obtained  as  part  of  the  solution  to  the  problem  and, 
indeed,  change  during  the  course  of  loading  and  unloading,  in  a  rather 
complex  manner.  Indeed,  in  our  model  the  tension  cracks  may  continue  to 
grow  at  a  state  of  unloading,  because  of  the  locking  of  the  preexisting 
flaw,  as  illustrated  in  the  examples  given  in  Sec.  4.  Furthermore,  we  have 
attempted  to  render  the  model  self-consistent  in  the  sense  that  the  stress 
intensity  factors  at  the  tips  of  the  tension  cracks  are  obtained  on  the 
basis  of  the  existing  gap  (created  due  to  the  sliding  of  the  preexisting 
flaw)  and  in  the  presence  of  farfield  loads,  and  these  are  then  made  to 
correspond  to  the  stress  intensity  factors  calculated  on  the  basis  of  the 
applied  loads  and  the  effective  total  force  transmitted  across  the 
preexisting  flaw.  In  the  model  of  Moss  and  Gupta,  the  stress  intensity 
factors  are  calculated  only  on  the  basis  of  the  applied  loads  and  the 
effective  force  across  the  preexisting  crack.  We  note  that  if  the 
preexisting  flaw  is  locked,  then  the  stress  intensity  factors  must  be 
computed  on  the  basis  of  the  existing  gaps  and  the  applied  loads,  and  these 
stress  intensity  factors  may  not  be  equal  to  the  ones  calculated  on  the 
basis  of  the  applied  loads  and  the  effective  force  across  the  preexisting 
flaw,  when  the  flaw  remains  dormant. 

3.  AVERAGE  QUANTITIES 

So  far,  the  problem  has  been  formulated  for  flaws  of  fixed  orientation 
4>  and  size  2c.  When  a  dilute  distribution  of  preexisting  flaws  is  involved, 
the  interaction  effects  may  be  neglected,  and  the  overall  strains  and  strain 
rates  may  be  estimated  by  a  simple  average  of  the  individual  contributions. 
To  this  end  let  N  -  N(^)  be  the  number  of  flaws  of  orientation  <f>  and  size  2c 
per  unit  area.  Hence,  f  -  c2N  is  given  as  a  function  of  4> .  The  average 


strain,  c,  then  may  be  computed  from 

£  -  -^-jy  £(rf.f(*))d*.  (3.1) 

where  c(^,f(^))  is  given,  in  our  analysis,  by  (2.8)  or  by  (2.27).  A  similar 
expression  can  be  written  for  the  average  strain  rate.  When  a  finite  number 
of  specific  orientations  of  preexisting  flaws  is  involved,  the  integral  in 
(3.1)  reduces  to  a  finite  sum, 

1  M 

T  --$r  l  «(*., f(*a>).  (3.2) 

o-l~ 

where  M  is  the  total  number  of  involved  orientations. 

In  the  formulation  presented  above  and  in  the  preceding  section,  no 
attempt  is  made  to  include  interaction  among  neighboring  flaws  and  their 
associated  cracks.  There  are  a  number  of  possible  ways  that  this 
interaction  may  be  included.  Perhaps  the  most  effective  way  would  be  to 
estimate  the  stress  intensity  factors  at  the  crack  tips  by  considering 
several  interacting  flaws  or  an  infinite  row  of  such  flaws,  in  the  manner 
discussed  by  Horii  and  Nemat-Nasser  (1983,  1985a, b,  1986).  This  poses  an 
extremely  complex  mathematical  problem  where,  although  solvable  as  shown  by 
Horii  and  Nemat-Nasser  (1985b),  the  resulting  equations  are  so  involved  that 
they  do  not  lend  themselves  to  simple  analytic  calculations  of  the  kind 
given  in  the  present  paper. 

An  alternative  approach  is  to  consider  an  averaging  technique  similar 
to  the  self-consistent  scheme;  see  Budiansky  (1965)  and  Hill  (1965).  This 
would  require  the  solution  of  the  problem  of  out-of-plane  crack  growth  from 
a  preexisting  flaw  embedded  in  a  linearly  elastic  but  anisotropic  material, 
which  currently  is  not  available.  Therefore,  any  improvement  on  the  simple 
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averaging  procedure  given  by  (3.1)  seems  to  require  either  extensive  and 
complicated  computations  beyond  the  scope  of  the  present  work  or  it  must  be 
based  on  additional  simplifying  assumptions,  which  is  not  our  intention  at 
this  time. 

In  the  following  section  we  shall  assume  a  dilute  distribution  of 
preexisting  flaws,  use  the  simple  averaging  given  by  (3.2),  and  present  a 
number  of  illustrative  examples. 

4.  RESULTS  AND  DISCUSSIONS 

In  this  section,  we  shall  illustrate  the  predictive  capability  of  the 
model  presented  in  the  preceding  section,  in  terms  of  several  numerical 
examples,  and  discuss  the  results. 

The  basic  parameters  of  the  model  are: 

1)  f  -  c2N  which  represents  the  density  and  the  size  of  the  flaws.  It  can 
be  given  any  suitable  distribution  in  terms  of  the  flaw  orientation  4> . 
In  our  examples  we  shall  assume  it  to  be  constant. 

2)  M0  which  is  the  number  of  flaw  orientations  between  0°  and  90°.  In  our 
examples  we  shall  assume  an  isotropic  initial  flaw  distribution,  and 
therefore  use  finite,  equally  spaced  orientations  between  0°  and  90°. 

3)  Kc  -  Kc/(r°7c)  which  is  the  dimensionless  critical  value  of  the  stress 
intensity  factor.  Since  typical  values  of  fracture  toughness,  Kc,  for 
rocks  are  of  the  order  of  105  -  10®  Nm'3/Z,  taking  a  typical  flaw  length 
to  be  of  the  order  of  10"*m,  and  the  cohesive  stress  to  be  of  the  order 
of  107  Nm"2,  we  estimate  Kc  to  be  of  the  order  of  1. 

4)  d0  which  represents  the  limiting  maximum  flaw  expansion  (dilatancy) 
normal  to  its  face;  see  Eq.  (2.5).  It  is  nondimensionalized  using  half 
the  flaw  length  c.  It  is  also  normalized  by  dividing  by  r°//j.  Hence, 


d0  -  1,  for  example,  represents  an  actual  maximum  possible  flaw  expansion 
of  0.5%  of  half  the  flaw  length  c  when  t°/h  -  5.0  x  10'3. 

5)  f  which  is  the  coefficient  of  the  exponential  saturation  of  the  flaw 

dilatancy  with  respect  to  the  flaw  slip.  Since  we  normalize  the  slip 
b  (and  the  dilatancy  d)  by  measuring  it  in  units  of  r3/V,  f  (and  17  in 
item  6  below)  is  measured  in  units  of  which  is  of  the  order  of 

102  -  103,  a  rather  large  dimensionless  number.  For  all  examples,  we  use 
f  -  1;  hence,  full  flaw  expansion  is  attained  at  flaw  slips  of  less  than 
1%  of  the  flaw  length. 

6)  t)  which  is  the  coefficient  of  the  exponential  decay  of  the  rate  of 
cohesive  resistance  of  the  preexisting  flaw. 

7)  Ji  which  is  the  coefficient  of  internal  friction  possibly  having  a  value 
between  0  and  1. 

8)  1/  which  is  the  Poisson  ratio  of  the  elastic  matrix.  In  actual  cases,  u 
should  be  interpreted  as  the  average  overall  Poisson  ratio  of  a  solid 
with  microdefects. 

For  the  examples  displayed  in  Figs.  3-6,  we  have  used  f  -  0.3,  Mg  -  10, 
v  -  0.25,  f  -  1.0,  and,  unless  otherwise  indicated  explicitly,  Kc  -  1.0, 
d0  -  0.05,  r)  -  0.05,  and  /I  -  0.30.  In  these  figures,  stresses  are 
normalized  by  r°  and  strains  by  (r3/^) . 

Figures  3a  and  b  show  typical  responses  in  uniaxial  compression.  In 
Fig.  3a  the  normalized  overall  volumetric  strain  and  in  Fig.  3b  the 

normalized  overall  axial  compressive  strain  are  plotted  on  the  horizontal 

axis  in  terms  of  the  normalized  compression  as  ordinate.  As  is  common  in 

the  rock  mechanics  literature,  compression  and  contraction  are  viewed  as 

positive.  Up  to  point  A,  the  overall  response  is  linearly  elastic,  because 
the  local  stresses  acting  on  the  preexisting  flaws  are  not  as  yet  large 


enough  to  activate  these  flaws.  At  point  A,  frictional  sliding  initiates  on 
some  flaws,  and  at  point  B,  cracks  nucleate  at  some  suitably  oriented  flaws. 
Upon  further  loading,  tension  cracks  begin  to  grow,  leading  to  substantial 
dilatancy  which  quickly  offsets  the  elastic  volumetric  contraction,  leading 
to  overall  volumetric  expansion  as  the  axial  compressive  load  is  further 
increased;  see  Fig.  3a. 

The  stress-strain  relation  shown  in  Fig.  3b  remains  monotonic  in 
uniaxial  loading,  and  closely  resembles  the  experimentally  observed  behavior 
of  rocks;  see,  e.g.,  Brace  and  Bombolakis  (1963).  In  the  presence  of 
lateral  tension,  however,  the  curve  would  terminate  at  a  point  when  a  most 
critical  tension  crack  begins  to  grow  in  the  axial  direction  in  an  unstable 
manner,  leading  to  axial  splitting;  see  Nemat-Nasser  and  Horii  (1982).  In 
the  presence  of  lateral  compression,  the  stress -strain  curve  remains 
monotonic  and  no  failure  is  predicted  by  this  model  because  the  interaction 
effects  among  neighboring  flaws  are  not  incorporated  in  the  present  model. 
As  has  been  shown  by  Horii  and  Nemat-Nasser  (1985a, 1986) ,  the  interaction 
among  a  row  of  suitably  oriented  preexisting  flaws  can  lead  to  the  formation 
of  a  fault,  as  tension  cracks  nucleate  at  these  flaws  and  grow  in  an 
unstable  manner.  The  analytical  computations  associated  with  such  a  model 
of  faulting  are  rather  complex  and  therefore  considerable  simplification  is 
necessary  before  such  a  failure  model  can  be  incorporated  into  a 
micromechanically  based  constitutive  relation  of  the  kind  developed  in  the 
preceding  section.  The  present  constitutive  model,  therefore,  is  not 
suitable  for  such  loading  regimes. 


At  point  C  in  Figs.  3a  and  b,  unloading  has  been  initiated.  The 
response  from  C  to  D  is  linearly  elastic.  However,  the  material  is  slightly 


more  compliant,  since  microcracking  has  led  to  certain 


stiffness 


degradation.  The  difference  between  the  elastic  moduli  at  initial  loading 
and  at  unloading  from  point  C  is  rather  small  in  the  axial  direction,  but  it 
is  considerable  for  lateral  tension  (but  not  for  lateral  compression).  This 
is  because  during  axial  compression,  cracks  grow  essentially  in  the  axial 


direction,  and  such  cracks  do  not  influence  much  the  axial  stiffness  of  the 
material,  but  they  do  have  considerable  weakening  effect  for  lateral 
tension. 


The  instantaneous  compliances  at  several  states  during  the  course  of  a 
monotonic  axial  compression  and  of  unloading  are  given  in  Table  1.  As  is 
seen,  strong  anisotropy  develops  due  to  slip  and  microcracking. 


From  point  C  to  point  D  in  Figs.  3a  and  b,  all  the  flaws  are  locked  and 
all  the  tension  cracks  remain  open.  Backsliding  of  some  preexisting  flaws 
initiates  at  point  D.  In  addition,  some  tension  cracks  grow  before  being 
dormant.  Considerable  volumetric  contraction  begins  to  occur  after  this 
state.  In  this  example,  no  crack  closure  occurs  during  the  entire  unloading 
process  which  leads  to  point  E  corresponding  to  some  permanent  total 
volumetric  expansion.  This  type  of  response  has  been  inferred  by  Scholz  and 
Kranz  (1974)  on  the  basis  of  experimental  observations. 


In  Figs.  4,  we  show  examples  of  loading-path-dependency.  Fig.  4a  shows 
the  stress -strain  curves  for  two  different  loading  paths  defined  in  Fig.  4b. 
On  path  I,  uniaxial  compression  is  applied  until  point  A,  producing  slip  and 
microcracks  at  some  preexisting  flaws.  Then,  with  axial  compression  held 
fixed,  the  lateral  compressive  stress  is  applied  up  to  point  C.  During  this 
regime,  no  further  slip  occurs  on  any  flaws  because  the  driving  shear 
stresses  on  the  flaws  actually  decrease,  but  not  enough  to  initiate 
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backsliding.  At  point  B,  closure  of  some  tension  cracks  begins. 


Path  II,  on  the  other  hand,  defines  a  proportional  loading  directly  to 
point  C' .  Unlike  for  the  path  I  loading,  for  path  II  no  slip  is  activated. 
The  response,  therefore,  is  linearly  elastic.  The  microstructures  of  the 
material  for  points  C  and  C'  are  quite  different. 

Our  model  includes  a  number  of  important  parameters,  and  we  now  proceed 
to  examine  their  effects  on  the  model's  predictions.  These  parameters  are 
listed  at  the  beginning  of  this  section.  Except  for  C.fj,  and  d0,  they  all 
have  rather  clear  physical  significance,  and  therefore  their  values  may  be 
estimated  on  the  basis  of  experimental  observations.  The  parameter  d0 
represents  maximum  expansion  (due  to  asperities  or  other  effects)  that  the 
flaw  itself  can  undergo  because  of  sliding.  The  parameters  f  and  r/ 
represent  the  rate  at  which  flaw  expansion  saturates  and  the  rate  at  which 
the  cohesive  resistance  of  the  flaw  decays,  both  measured  per  unit  rate  of 
sliding.  These  parameters  have  a  rather  significant  influence  on  the 
overall  volumetric  response  predicted  by  the  model. 

Figures  5  and  6  show  the  effect  of  ij  and  of  d0  in  uniaxial  compression, 
followed  by  unloading.  As  pointed  out  before,  when  (r°/p)  -  5.0  x  10~3, 
d0  -  1  represents  an  actual  maximum  limiting  flaw  expansion  normal  to  its 
surface  of  only  0.5%  of  half  the  total  flaw  length,  and  hence  it  is  small. 
Nevertheless,  it  has  a  rather  significant  influence  on  the  response 
predicted  by  the  model. 

The  influence  of  internal  friction  Ji  and  the  normalized  fracture 
toughness  Kc  on  the  dilatancy  in  uniaxial  compression  followed  by  unloading 
is  depicted  in  Figs.  7  and  8,  respectively.  As  one  would  expect,  a  greater 


dilatancy  results  at  lower  frictional  resistance  of  the  flaws  and  lower 
fracture  toughness  of  the  matrix  material.  It  should  be  noted  that  there  is 
a  considerable  difference  between  the  effects  of  the  frictional  parameter  Ji 
and  the  normalized  fracture  toughness  Kc  on  the  overall  response:  whereas 
the  friction  coefficient  jl  affects  the  flaw  response  during  the  entire 
loading  and  unloading  regime,  Kc's  influence  emerges  only  with  tension  crack 
initiation.  Note  also  that  normalized  Kc  involves  the  flaw  size  2c  and  the 
cohesive  stress  r°,  as  well  as  the  actual  fracture  toughness  of  the 
surrounding  material.  Therefore,  while  the  order  of  magnitude  of  Kc  can  be 
estimated  on  the  basis  of  experimental  observations,  its  actual  value 
involves  quantities  that  are  difficult  to  measure.  For  example,  the 
fracture  toughness  should  actually  reflect  the  behavior  of  aggregates  of 
crystals  involving  a  variety  of  preexisting  defects  in  a  rock  or  ceramic 
sample . 

As  a  final  illustration,  we  seek  to  fit  the  model  to  experimental 
results  reported  by  Zoback  and  Byerlee  (1975) .  Since  actual  rocks  contain 
preexisting  open  cracks  and  pores  which  tend  to  close  upon  application  of 
loads,  the  initial  stress-strain  relations  often  reported  are  considerably 
softer  than  the  subsequent  response.  We  have  not  included  this  in  our 
calculations,  although  it  can  easily  be  done,  as  detailed  by  Horii  and 
Nemat-Nasser  (1982).  Therefore,  we  do  not  seek  to  match  the  initial  portion 
of  the  stress-strain  curves,  but  only  compare  their  overall  shapes.  To 
offset  the  initial  soft  response,  we  simply  shift  the  model's  prediction,  as 
shown  in  Figs.  9a  and  b. 

Since  the  model  is  a  two-dimensional  one,  whereas  the  experimental 
results  are  for  uniaxial  compression  of  a  cylindrical  sample,  a  certain 


amount  of  adjustment  in  the  model  parameters  is  required.  The  reported  flaw 
size  2c  is  of  the  order  4  ~  8  x  10'4,  and  taking  f  -  N(c/V)3  -  0.2,  for  the 
two-dimensional  case  we  choose  f  -  (0.2)2/3  -  0.34.  The  Poisson  ratio  and 
the  shear  modulus  are  estimated  to  be  v  -  0.35  and  fi  —  2.1  x  104  MNm"2, 
based  on  the  data  reported  by  Zoback  and  Byerlee  (1975) .  In  addition,  we 
take  Kc  -  0.7  MNnT3/2,  r°  -  21  MNm'2,  and  set  2c  -  5.0  x  10"4m  (therefore, 

Kc  -  2.1).  The  coefficient  of  friction  is  taken  as  a*  -  0.1,  which  is  rather 
small,  but  may  be  justified  on  the  grounds  that  many  of  the  flaws  seem  to 
have  been  initially  open  and  therefore  their  frictional  resistance  ought  to 
be  small.  Considering  the  simplicity  of  the  model,  the  fit  shown  in  Figs. 
9a  and  b  is  rather  good.  Nevertheless,  in  view  of  the  number  of  free 
parameters  involved,  the  comparison  should  be  regarded  as  an  illustration 
and  not  necessarily  a  direct  support  of  the  viablity  of  the  model. 

In  summary,  it  may  be  concluded  that  the  micromechanical  model 
consisting  of  preexisting  flaws,  with  frictional  and  cohesive  resistance  to 
sliding,  and  capable  of  nucleating  tension  cracks,  seems  to  capture  a  number 
of  important  features  observed  in  the  response  of  brittle  materials  with 
microdefects.  The  model  presented  here  does  not  include  the  important 
interaction  effects  among  neighboring  flaws.  Judging  from  results  reported 
by  Horii  and  Nemat-Nasser  (1986)  on  the  failure  in  compression  of  materials 
of  this  kind  under  moderate  and  large  confining  pressures,  we  feel  that 
constitutive  modeling  based  on  the  mechanisms  considered  by  these  authors 
should  prove  to  be  effective. 


Figure  Captions 


Figure  1:  An  isolated  preexisting  flaw  PP'  with  tension  cracks  P'Q'  and 
PQ. 

Figure  2:  (a)  An  equivalent  isolated  crack  with  a  pair  of  colinear 
concentrated  forces  applied  at  its  center.  (b)  An  equivalent 
isolated  crack  with  concentrated  gap  at  its  end. 

Figure  3:  (a)  Volumetric  strain  versus  differential  stress  curve  for  a 

loading  and  unloading  loop.  (b)  Axial  strain  versus  stress 
curve  for  a  loading  and  unloading  loop;  f  -  0.3,  M0  -  10, 
i/  -  0.25,  £  -  0.3,  Kc  -  1.0,  dQ  -  0.05,  and  r,  -  0.05. 

Figure  4:  (a)  Axial  strain  versus  differential  stress  for  different 
loading  paths;  f  -  0.3,  Mg  -  10,  i/  -  0.25,  m  —  0.3,  d0-0, 
Kc  —  1.0,  and  r)  —  0.05.  (b)  Loading  paths. 

Figure  5:  Volumetric  strain  versus  differential  stress  curves  for  various 
rj;  f  -  0.3,  M0  -  10 ,  i/  -  0.25,  p  -  0.3,  d0-0,  and  Ke  -  1.0. 

Figure  6:  Volumetric  strain  versus  differential  stress  curves  for  various 
d0 ;  f  -  0.3,  M0  -  10,  v  -  0.25,  m  -  0.3,  Kc  -  1.0,  and  rj  -  0.05. 

Figure  7;  Volumetric  strain  versus  differential  stress  curves  for  various 
MI  f  -  0.3,  M0  -  10,  „  -  0.25,  Kc  -  1.0,  d0  -  0.05,  and 

t)  -  0.05. 

Figure  8:  Volumetric  strain  versus  differential  stress  curves  for  various 
Kc ;  f  -  0.3,  M0  -  10,  -  0.25,  Ji  -  0.3,  d0  -  0.05,  and 

rt  -  0.05. 
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APPENDIX  A:  Strains  Due  to  Crack  Opening 


The  strain  due  to  the  crack  opening  displacement  is  given  by 

£ij  -y  y[  t  ni  t  uj  ]  +  nj  [  ui  ] )  ds  ,  (A  •  1) 

J  ~C 

where  V  is  the  volume,  N  is  the  number  of  cracks  in  this  volume,  n  is  a  unit 
vector  normal  to  a  crack  surface,  and  [uj  is  the  i-th  component  of  the 
crack  opening  displacement.  From  this,  the  inelastic  strain  due  to  crack 
slip  is  computed  as  follows.  In  a  local  coordinate  system  with  the  crack  on 
the  Xj-axis, 


£12  “  e2i  "  Y  I  tu2ldx  -  Ybc  "  fb’ 

J  0 

£n  “  £22  ”  0|  b  —  b/c . 


(A. 2) 


Transformation  to  the  global  coordinate  system  is  made  with  the  following 
matrix: 


T 


cos  d  sin  d 
-sin  d  cos  d  ’ 


(A. 3) 


where  d  is  the  angle  between  the  local  and  global  x^axis.  By  use  of  this 
matrix,  we  have 


£b  " 


-sin  2^  cos  2d 
cos  2 d  sin  2d  ’ 


(A. 4) 


in  the  global  coordinate  system. 


For  the  Mode  I  crack  displacement,  we  have 


-li 


2  x 


•I 


[u,]dx  -  2fd, 


(A. 5) 


in  the  local  coordinates,  and 


£d  “  2fd 


(A. 6) 


cos2^  cos  ^sin  <t> 
cos  ^sin  4>  sin2^ 

in  che  global  coordinates . 
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APPENDIX  B:  Strains  Due  to  a  Dislocated  Crack 


First,  we  examine  the  Mode  I  crack  displacement.  Consider  an  array  of 
edge  dislocations  continuously  distributed  along  -a<x2<a,  x1-0;  see 
Fig.  B.l.  The  corresponding  normal  stress  is  (see,  e.g.,  Muskhelishvili , 
1953). 

i  u  r*  Bi(o 

■’n<**.°>  -  -S-T^T  <B1> 

and  the  displacement  in  the  Xj^-direction  on  the  x2-axis  is  given  by 

rx2 

u1(0,x2)  -  g(x2)  -  -  Bx(Od£,  (B.2) 

-  oo 

where  B1(x2)  is  the  dislocation  density.  Since  crack  surfaces  are 
traction- free ,  with  the  farfield  stresses  applied,  the  following  integral 
equation  defines  the  dislocation  density: 

•  Bx<0 

7 - —  +  p(x2)  -  0,  -a  <  x2  <  a,  (B.3) 

.  S  *  x2 

where  p(x2)  is  the  traction  on  the  x2-axis  due  to  farfield  stress. 

A  solution  of  (B.3)  singular  at  x2  -  ±a  is  given  by  (e.g., 
Muskhelishvili,  1953) 

/.[''f-  £  p(°d(  +  Cl’  <B'4) 

where  C  is  a  constant  yet  to  be  determined  by  another  boundary  condition. 
In  this  case,  the  farfield  stress  is  zero,  but  the  crack  is  dislocated  at 
x2  -  a, 

;  P(x2)  -  0,  [uj (a)  -  gj.  (B.  5) 

| 

From  (B.2),  (B.4),  and  (B.5),  we  have 

\ 

r 


2*  1  -  V 


Note  that  the  condition  [ ut ]  -  0  for  -a  >  x2  is  automatically  satisfied  by 
this  array  of  dislocations.  Substituting  (B.6)  into  (B.l),  we  obtain  the 
normal  stresses  on  xx  -  0, 

gi  u  sgn(x2) 

*u<x2.0)  “  -!7TrV  J4  -  a2  H (  1  x2 1  -  a),  (B.7) 

where  H(x)  is  the  Heaviside  step  function.  The  stress  intensity  factor  at 
x2  *  -a  now  is 


_  H  1 

1  2  1-i/  Jrta. 


(B .  8) 


From  (B.2),(B.5),  and  (A.l),  we  compute  the  strain  components  in  the  local 
coordinate  system  as  follows: 


£n  ”  6ia. 

£12  "  £21  ”  £22  “  0- 

Note  that  both  gx  and  a  are  dimensionless;  see  Appendix  A. 

Now  consider  the  Mode  II  crack  displacement.  If  the 
at  x2  *  a  is  denoted  by  g2,  we  obtain 


1 


£12  ”  £21  “  2  ®2a' 


(B .  9) 


concentrated  gap 


(B . 10) 


In  our  particular  problem, 
g:  -  b  sin  6  +  d  cos  6 , 


g,  -  b  cos  9  -  d  sin  9 , 


(B . 11) 


Thus,  considering  that  we  have  a  pair  of  dislocated  cracks,  the 
components  in  the  local  coordinate  system  are  given  by 


2sin 

9 

cos  9 

*  di 
+  f  — r— 

2cos  9 

-sin  9 

cos 

9 

0 

2 

-sin  9 

0 

The  transformation  to  the  global  coordinate  system  is  made  with 


cos  (0  +  <f>)  sin  ( 9  + 

-sin  (9  +  cos  (9  +  <(>)  ’ 


and  we  obtain 


£c  -  fbiy 


■2sin  ^cos (9  +  <f>) 
cos  (9  +  2  <f>) 


cos  (9  +  2  <f>) 
2cos  <f>sin(9  +  <j>) 


+  fdiy 


2cos  <f>cos(9  +  <f> ) 
sin(0  +  2  <f>) 


strain 


(B.12) 


(B. 13) 


(B. 14) 


in  the  global  coordinate  system. 


sin(0  +  2<t> ) 
2sin  $sin(0  +  <j>) 
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APPENDIX  C:  Strains  Due  to  Crack  Displacement  Caused  by  Applied  Stresses 

Under  the  farfield  stresses  T  and  S,  the  crack  opening  displacements 
are  given  by 


[uj  -  2T  1  ’  —  J{i/ 2)2-x|, 


[u2]  -  2S 


M 

1  -  1/ 


(C.l) 


y(i/2)2-x2. 


In  our  particular  problem,  T  -  a:a  and  S  -  o:0/2.  Therefore,  after 
integrating,  the  strains  in  the  local  coordinate  system  are  expressed  as 


£a 


1  - 
8M 


:  a  o\fi 
a:) 9  ~0~ 


(C .  2) 


Using  (B.ll),  we  have  the  expression  in  the  global  coordinate  system, 
£.  “  f  ~s/  **2[4(q:a)g  +  (a:g)g]. 


(C .  3) 


APPENDIX  D:  Stress  Intensity  Factor  of  a  Dislocated  Crack 


Consider  the  case  where  the  crack  opening  displacement  normal  to  the 
flaw  is  zero.  By  use  of  (B.8),  the  Mode  I  stress  intensity  factor  of  a 
dislocated  crack  of  length  2  is  obtained  to  be 


K 


I 


b 


+  a:a  J 


irZ 
2  • 


(D.l) 


Since  this  expression  is  singular  when  2  -  0,  we  replace  it  by  the  following 
non- singular  form 


V  -  1  -  : - - -  b  +  a:aj 

1  '  u  J2x(2  +  r*> 


7T  j! 


(D.2) 


and  seek  to  estimate  i**  such  that  at  2-0*,  Kf  -  Kr.  When  2-0,  the 
average  slip  of  the  preexisting  flaw  equals  the  average  Mode  II  crack 
displacement,  and  is  given  by 

"  “j^T  |  2<2‘-E2  *  Te>  ~ '  J*2  '  *2  dx  ~  Y~~jt  ^  *c(2:E2 

"C 

On  the  other  hand,  when  2—0,  (2.9a)  reduces  to 


rc)  sin  9, 


which  is  maximum  at  6  -  x/2 . 
is  determined  to  be 


32 


2*. 


From  (D.2), 


(D.4) 

(D. 3) ,  and  (D.4)  with  6  -  n/2 ,  2 ** 


(D.5) 


We  simply  generalize  this  result  when  dilatancy  d  is  included.  The  final 
form  of  Kj'  then  is 


V  - 


1  -  V 


J2x(2  +  2 **) 


(b  sin  9  +  d  cos  9)  +  a:aj 


» 2 


(D.6) 
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APPENDIX  E:  Some  Basic  Equations  in  Loading 


•‘v'wlv'v!? 


'■I  n 


-77  —  -  —  -  (a: p,  -  r.)  sm 

dl  2  {*(i  +  i*)}3/2  C2 


w  I  z  / 
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77T 
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(sin  0  +  f(d0  -  d)  cos  0), 
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_a_  SKj' 
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2  1  -  ^  {2 w(i  +  i**)}3'2 


(b  cos  0  -  d  sin  9)  -  a 
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APPENDIX  G:  Stress  Intensity  Factor  in  Unloading 


Equation  (2.25)  with  Kt  defined  by  (2.24a), 
may  be  explained  in  the  following  way. 

At  the  end  of  a  loading  regime  when 
the  applied  farfield  stresses  are  £7nax, 
the  residual  stress  on  the  preexisting  flaw 
is  given  by 

rm  “  £max  •  P  2  *  rc  (G.l) 

Here,  rm  is  the  maximum  shear  stress  experienced  by  the  flaw  during  this 
loading  regime.  Backsliding  is  considered  to  occur  when  this  residual 
stress  overcomes  the  frictional  and  cohesive  resistance  of  the  flaw,  i.e., 
when 


K{  -  Kj.  (G . 4) 

The  essential  difference  between  (2.24a)  and  (2.10a)  is  due  to  the  t  rms 
associated  with  the  farfield  applied  stress,  i.e.,  o:aJ(ir2/2)  and  o-.aj(iri). 
This  difference  stems  from  the  fact  that  the  computation  of  the  stress 
intensity  factor  in  the  presence  of  a  gap  is  for  a  crack  of  length  i, 
whereas  that  associated  with  the  force  F  is  for  a  crack  of  length  2i. 
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Because  tension  cracks  grow  in  the  the  maximum  compressive  stress  direction, 
this  difference  is  quite  small.  Even  for  extreme  cases,  e.g.,  just  after 
crack  nucleation,  i  is  small  and  hence  this  term  is  very  small.  We  may 
therefore  use  the  expression  (2.25)  and  (2.26)  with  Kx  defined  by  (2.24a), 
in  order  to  ensure  a  consistent  formulation  without  loss  of  accuracy. 


APPENDIX  H:  Partially  Closed  Dislocated  Crack 


We  consider  the  case  shown  in  Fig.  B.l.  In  this  case,  the  integral 
equation  corresponding  to  Eq.  (B.3)  becomes 

*  Bx<0 


J. 


t  -  x 


d£  -  2 *  1-~  T. 


(H.l) 


If  we  solve  this  equation  under  the  condition  that  ctu  is  bounded  at  x2  -  -a 
and  unbounded  at  x2  -  a,  we  obtain  (e.g.,  Muskhelishvili ,  1953) 

i-i/  y®  x2 


Bx(x)  -  2T 


7 


a  -  x. 


(H.  2) 


After  integration  of  (B.2),  the  crack  opening  displacement  is  given  by 


[UxHXj)  -  - 


The  condition  at  x2  -  a  yields 


—  2T  [  a  s  in' 1  —  +  a  —  -  Ja 2  -  x|  ] 


(H.  3) 


a  -  - 


jt _  _§i 


2*  1  -  u  T 


(H.4) 


Apparently,  this  result  is  reasonable  only  when  T  <  0.  In  our  particular 
problem, 


gj  -  b  sin  8  +  d  cos  8,  T  -  o:a.  (H.5) 

Thus,  integration  of  the  crack  opening  displacement  leads  to  the  following 
expression  for  the  strain  component  cu  in  the  local  coordinate  system: 

en  -  -f  -  ,  1  v  (b  sin  8  +  d  cos  8)2  .  (H.6) 

11  2n  1  -  v  (o:a) 

The  strain  component  e12  requires  no  modification  when  we  assume  no 
friction  on  tension  cracks.  Considering  the  contribution  by  the  farfield 
stresses  as  well  as  the  dislocation  array,  we  obtain 


Is 

*  S' 

1 

v!v 


,:k 

§ 

3$ 


I 

v.V 

vS 

$$ 

'M 


£«  -  f  T  (b  cos 


d  sin  *)  +  f-~'---*i«(g;g) 


(H.  7) 


in  the  local  coordinate  system.  The  transformation  to  the  global  coordinate 
system  is  done  by  Eq.  (B.13).  In  the  global  coordinate  system,  we  hence 


-  f  ^ r  (b  sin  6  +  d  cos  0)2a 

2k  1  -  v  ~ 


+  f  -r-  (b  cos  8  -  d  sin  8)p 


(H.  8) 
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Table  1:  Instantaneous  compliance,  at  various  loading  and  unloading  states. 


an/r° 

cu  /(r2/^) 

i 

^D1122 

4D2222 

PD2211 

4^1212 

0* 

0 

0.375 

;  -0.125  1 

0.375  : 

-0.125 

0.25 

3.68 

1.42 

0.421 

-0.212  ! 

0.462  ; 

-0.171 

0.257 

7.46 

3.21 

0.483 

-0.652  i 

3.775  1 

-0.510 

0.331 

13.  r* 

5.99 

0.391 

-0.391  ! 

27.6 

-1.18 

1  2.368 

5.28* 

2.95 

0.472 

-1.03 

14.2 

-0.648  ! 

2.657 

*  Initial  elastic  loading 
**  Upon  unloading 
f  Unloading  with  backsliding 
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